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Abstract. In this paper, we initiate the rigorous mathematical study of the problem of 
impulsive gravitational spacetime waves. We construct such spacetimes as solutions to the 
characteristic initial value problem of the Einstein vacuum equations with a data curvature 
delta singularity. We show that in the resulting spacetime, the delta singularity propagates 
along a characteristic hypersurface, while away from that hypersurface the spacetime re- 
mains smooth. Unlike the known explicit examples of impulsive gravitational spacetimes, 
this work in particular provides the first construction of an impulsive gravitational wave 
of compact extent and does not require any symmetry assumptions. The arguments in the 
present paper also extend to the problem of existence and uniqueness of solutions to a larger 
class of non-regular characteristic data. 
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1. Introduction 

In this paper, we initiate the study of the characteristic initial value problem for impulsive 
gravitational spacetimes in general relativity. These were considered to be solutions of the 
vacuum Einstein equations 

R^v = 

with a delta singularity in the Riemann curvature tensor supported on a null hypersurface. 
Their historical origin can be traced back to the cylindrical waves of Einstein- Rosen [15] , the 
plane waves of Brinkmann [7J and the explicit impulsive gravitational spacetimes of Penrose 



Impulsive gravitational waves have often been studied within the class of plane fronted 
gravitational waves. This class of explicit solutions to the vacuum Einstein equations has 
been first studied by Brinkmann [7] and interest in them has been revived in later decades by 
the work of Bondi-Pirani- Robinson [6 J . They were later classified geometrically by Jordon- 
Ehlers-Kundt [20]. Among this class are the pp-waves (plane fronted waves with parallel 
rays) that were discovered by Brinkmann [7J, for which the metric takes the form 

g = Idudr + H(u, X, Y)du 2 + dX 2 + dY 2 , 

and the Einstein vacuum equations imply that 

d 2 H d 2 H 



Q X 2 + dY 2 °' W 

These include the special case of sandwich waves, where H is compactly supported in u. 
Since ([I]) is linear, it follows that pp-waves enjoy a principle of linear superposition. 

pp-w&ves have a plane symmetry and originally impulsive gravitational waves have been 
thought of as a limiting case of the pp-wave with the function H admitting a delta singularity 
in the variable u. Precisely, explicit impulsive gravitational spacetimes were discovered and 
studied by Penrose [35], who gave the metric in the following double null coordinate form: 

g = -2dudu + (1 - uQ(u))dx 2 + (1 + u<d(u))dy 2 , 

where O is the Heaviside step function. In the Brinkmann coordinate system, the metric has 
the pp-wave form and an obvious delta singularity: 



g = -2dudr - 5(u)(X 2 - Y 2 )du 2 + dX 2 + dY 2 , 
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where 5(u) is the Dirac delta. Despite the presence of the delta singularity for the metric 
in Brinkmann coordinates, the corresponding spacetime is Lipschitz and it is only the Rie- 
mann curvature tensor (specifically, the only non-trivial a component of it) that has a delta 
function supported on the plane null hypersurface u — 0. This spacetime turns out to pos- 
sess remarkable global geometric properties [31]. In particular, it exhibits a strong focusing 
property such that the whole null cone emanating from a point in the past of u = refocuses 
to a single line in the future of u = (See Figure 1). In [34J, this property forms the basis 
of Penrose's argument that global hyperbolicity fails in this spacetime. Penrose's impulsive 
spacetime is plane symmetric, non-asymptotically flat and the delta curvature singularity 
is supported on the 3-dimensional infinite plane {(u, u,x, y) : u = 0}. It has been long 
debated whether the strong focussing property and the resulting lack of global hyperbolicity 
is directly tied to the infinite extent of the impulsive gravitational wave (See Yurtsever [43]). 




Figure 1. Focusing in Penrose's impulsive gravitational spacetime 



While the interest in impulsive gravitational spacetimes had been high due to the avail- 
ability of explicit solutions, their global geometric properties, superposition properties and 
their limiting relation to general spacetimes (see Penrose [35] , Aichelburg-Sexl [TJ), the 
first study of general spacetimes satisfying the Einstein equations and admitting possible 
three dimensional delta singularities was undertaken by Taub [12], who derived a system of 
consistency relations linking the metric, curvature tensor and the geometry of the singular 
hypersurface. Impulsive gravitational spacetimes also arise as high-frequency limiting space- 
times considered by Choquet-Bruhat [8]. We refer the readers to [TE] . [IT] . [3J, [5] and the 
references therein for a more detailed exposition on the physics literature. 

In this paper, we begin the study of impulsive gravitational spacetimes viewed in the 
context of an (characteristic) initial value problem. We consider the data, prescribed on an 
outgoing null hypersurface u = and assume that the a component of the curvature has 
a delta singularity supported on a two dimensional surface S 0>u , which can be thought of 
as the intersection of the hypersurfaces u = and u = u s . Observe that in the Penrose's 
explicit impulsive solution, the a curvature component has precisely this type of behavior 
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when restricted to u — 0. Unlike the explicit impulsive spacetimes, which have only been 
constructed in plane symmetry and thus have infinite spatial extent, we will consider the 
case of So t u s with compact topology, more precisely a sphere. The data on the incoming null 
hypersurface is prescribed to be smooth but otherwise without any smallness assumptions. 

With this data, we show that a unique spacetime satisfying the vacuum Einstein equations 
in an appropriate sense can be constructed locally. Moreover, the delta singularity propagates 
along a null hypersurface emanating from the initial singularity on So,u s an d the spacetime 
is smooth away from this null hypersurface (See Figure 2). 




So,o 



Figure 2. Propagation of Singularity 

Our main result on the propagation of an impulsive gravitational wave is described by 

Theorem 1. Suppose the following hold for the initial data set: 

• The data on H_ is smooth. 

• The data on H is smooth except across a two sphere So jU , where the traceless part 
of the second fundamental form x has a jump discontinuity. 

Then 

(a) Given such initial data and e sufficiently small, there exists a unique spacetime (Ai,g) 
endowed with a double null foliation u, u and solves the characteristic initial value prob- 
lem for the vacuum Einstein equations in the region < u < «*, < u < whenever 
u*,u* < e. 

(b) Define H_ u to be the incoming null hypersurface emanating from So M . Then the cur- 
vature components qiab = R(^A, 64, ee, 64) are measures with a singular atom supported 
on the hypersurface H_ u . All other components of the curvature tensor can be defined in 
L°° . Moreover, the solution is smooth away from H_ u . 

For an impulsive gravitational wave, standard local existence shows that the spacetime is 
smooth in < u < u s . The problem of constructing an impulsive gravitational wave lies in 
making sense of the solution in the whole region with u s < u < and showing that the 
singularity propagates along the characteristic hypersurface H_ u and that the spacetime is 
smooth after the impulse. 

Remark 1 (Larger class of initial data). The proof introduces a new type of energy estimates 
for the Einstein vacuum equations which allows some components of the Riemann curvature 
tensor not to be in L 2 . This discovery allows us to consider the problem of local existence and 
uniqueness for a larger class of non-regular data. This includes data with Riemann curvature 
tensor that can only be understood as a conormal distribution. 

Remark 2 (Uniqueness). Penrose's construction of explicit impulsive gravitational space- 
times is based on the gluing approach, which has later been used to generate other explicit 
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solutions of the Einstein equations. In the general case of the characteristic initial value 
problem for an impulsive gravitational wave, an appropriate adaptation of the gluing phi- 
losophy would allow us to construct weak solutions with undetermined uniqueness. In this 
paper, however, we construct strong solutions in the sense that that uniqueness can also 
be established. To prove uniqueness, we establish a priori estimates for a larger class of ad- 
missible initial data. This allows us to show that the solution is unique among all C° limits 
of smooth solutions to the vacuum Einstein equations. 

Remark 3 (Propagation of singularity). Our theorem gives a precise description of the 
propagation of the initial singularity. Such problems have been extensively studied for semi- 
linear equations (see for example [I], [S^. Our result can be formally compared to the works 
of Majda on the propagation of shocks [2S], [SO] for systems of conservation laws and the 
subsequent [2], [52], [S3]? [II]- In these works, a short time existence, uniqueness and regu- 
larity result was established for initial data with a jump discontinuity across a surface and 
a precise description of the propagation of the singularity was also given. We establish an 
analogous result of propagation of singularity for a nonlinear system of quasilinear hyperbolic 
equations. However, contrary to [29], [30] the singularity that is considered in the present 
paper is not a shock, as it propagates along the characteristics. Moreover, unlike in [29] , [30J ; 
[2j ; [32j, |33j ; [H], where the problem is reformulated as an initial-boundary value problem 
and uniqueness is known only within the class of piecewise smooth solutions, our solution is 
also unique among limits of smooth solutions. In order to achieve this, the special structure 
of the Einstein equations in the double null foliation gauge has been heavily exploited. 

The construction of spacetimes from non-regular characteristic initial data consistent with 
that of an impulsive gravitational wave had been known only under symmetry assumptions 
[9], [26], [27]. The work of Christodoulou, who solved the characteristic initial value problem 
for data with bounded variation for the spherically symmetric Einstein-scalar field system 
[D], can be thought of as a first result in that direction. In particular, in this work, the 
second derivatives of the scalar field, which formally is analogous to the curvature, is allowed 
to have a delta singularity. The study of data with bounded variation turned out to have 
important consequences in the global structure of spacetimes and the resolution of the cosmic 
censorship conjecture for the spherically symmetric Einstein-scalar field system [TO], jllj . 
The construction of distributional solutions for the vacuum Einstein equations that include 
plane impulsive gravitational wave was carried out in [26] , [27] . The present paper is the first 
work that provides a consistent study of the initial value problem of impulsive gravitational 
spacetimes, including their existence, uniqueness and propagation of singularity /regularity. 

Going beyond spacetimes which represent a single impulsive gravitational wave, colliding 
impulsive gravitational waves had been studied by Khan-Penrose [21] and Szekeres [U]. 
In these explicit solutions, the spacetimes possess two null hypersurfaces with curvature 
delta singularity with a transverse intersection, representing the nonlinear interaction of two 
impulsive gravitational waves. The study of the characteristic initial value problem for the 
colliding impulsive gravitational waves will be carried out by the authors in a subsequent 
paper. 

1.1. First Version of the Theorem. Our general approach is based on energy estimates 
and transport equations in the double null foliation gauge. This general approach in the 
double null foliation gauge has been carried out in [22], [T2] and |24J. 
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The spacetime in question will be foliated by families of outgoing and incoming null 
hypersurfaces H u and H_ u respectively Their intersection is assumed to be a 2-sphere denoted 
by S Ui u- Define a null frame {e±, e 2 , e 3 , e 4 }, where e 3 and e 4 are null, as indicated in Figure 
3, and e±, e 2 are tangent to the two spheres S UM . e 4 is tangent to H u and e 3 is tangent to 




Figure 3. The Basic Setup and the Null Frame 

Decompose the Riemann curvature tensor with respect to this frame: 

a AB = R{e A) e 4 , e B , e 4 ), a AB = R(e A , e 3 , e B , e 3 ), 

Pa = ^R(ca, e 4 , e 3 , e 4 ), §_ A = ^R(e A , e 3 , e 3 , e 4 ), 

p = |-R(e 4 , e 3 , e 4 , e 3 ), a = \ * R ( e ^ e 3> e 4, e 3 ) 

In the context of impulsive gravitational spacetimes, the a component can only be understood 
as a measure. 

Define also the following Ricci coefficients with respect to the null frame: 

Xab = g(D A e 4 , e B ), x AB = g(D A e 3 , e B ), 

Va = -^g(D 3 e A , e 4 ), r[ A = -^g(D 4 e A , e 3 ) 
w = ~g{D A e 3 ,e^), u = -^g(D 3 e 4 ,e 3 ), 
(a = ^g(D A e 4 , e 3 ) 

Define x ( r esp. x) to be the traceless part of x (resp. x)- For the problem of the 
propagation of impulsive gravitational waves, we prescribe initial data on Hq such that x 
has a jump discontinuity across So jU but smooth otherwise. On H_ , we prescribe the initial 
data to be smooth but without any smallness assumptions. 

As mentioned before, we will prove local existence and uniqueness for a class of data more 
general than that for the impulsive gravitational wave. More precisely, we require that x 
and its angular derivatives are merely bounded. However, we do not require the derivative 
of x i n the e 4 direction even to be defined. 

Theorem 2. Suppose the characteristic initial data are smooth on if . On H Q , the charac- 
teristic initial data are determined by x and ( A . Let 6 be the coordinates on the two sphere 
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foliating H . Suppose, in every coordinate patch 

j<3 fc<l i<3 — 

Then there exists e sufficiently small such that the unique solution to the vacuum Einstein 
equations (Ai,g) endowed with a double null foliation u, u exists in < u < e, < u < e. 
Associated to the spacetime a coordinate system (u,u,9 1 ,9 2 ) exists, relative to which the 
spacetime is in particular Lipschitz and retains higher regularity in the angular directions. 

Remark 4. We will also explicitly construct a class of initial data that satisfy the assump- 
tions of Theorem [7] in Section [3| This construction of initial data will indicate an easy 
modification leading to an even more general class satisfying the assumptions of Theorem^ 

This Theorem in particular implies Theorem [jja). An additional argument, based on the 
estimates derived in the proof of Theorem [2j will be carried out in Section 5.10 to prove 



part (b) of Theorem [TJ The statement of Theorem [2] will be made precise in Section 2.5 In 
particular, the smoothness assumptions of Theorem [2] can be weakened and we will define 
in what sense the vacuum Einstein equations are satisfied. 

Remark 5 (Comparison to local existence results). Without symmetry assumptions, all 
known proofs of existence of spacetimes satisfying the Einstein equations are based on L 2 - 
type estimates for the curvature tensor and its derivatives or the metric components and their 
derivatives. Even with the recent resolution of the L 2 curvature conjecture by Klainerman- 
Rodnianski-Szeftel [25], [37], [38], [39], [10], the Riemann curvature tensor has to be at least 
in L 2 . (The classical local existence result of Hughes-Kato-Marsden [19] requires the metric 
to be in H s , with s>%.) 

In this paper, we prove local existence and uniqueness under the assumption that the space- 
time is merely Lipschitz, which in terms of differentiability is even one derivative weaker than 
the L 2 curvature conjecture. Of course the Lipschitz assumption refers to the worst possible 
behavior observed in our data and our result heavily relies on the structure of the Einstein 
equations which allows us to efficiently exploit the better behavior of the other components. 

1.2. Strategy of the Proof. For an impulsive gravitational wave, the curvature tensor can 
only be defined as a measure and is not in L 2 . This is one of the main challenges of this 
work. Let \I> denote the curvature components and ip denote the Ricci coefficients. The key 
new observation of this paper is that the L 2 -type energy estimates for the components of the 
Riemann curvature tensor 



H u Jh^ J h Jh Jo Jo Js uI%i 
coupled together with the null transport equations for the Ricci coefficients 

V 3 ip = * + ipip, V 4 V> = * + ipip 
can be renormalized and closed avoiding the L 2 -non-integrable components of curvature. 

1.2.1. Renormalized Energy Estimates. The difficulty in carrying out the above argument is 
that in our setting, not all curvature components are defined in L 2 . Even at the level of initial 
data, a is only defined as a measure. Therefore, we need to prove L 2 energy estimates without 
involving a. To this end, we need to introduce and estimate the renormalized curvature 
components. The idea of renormalizing the curvature components has been introduced in 
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j. Unlike this work, there the renormalization was used to obtain estimates for the Ricci 
coefficients. 

The classical way to derive energy estimates for the Einstein equations is via the Bel- 
Robinson tensor. In view of our renormalization, we avoid the use of the Bel-Robinson 
tensor and instead prove the energy estimates directly from the Bianchi identities (|8]). We 
note that the derivation of the energy estimates directly from the Bianchi equations without 
using the Bel Robinson tensor has also appeared in the work of Holzegel [IE] in a different 
setting. The challenge and motivation of this method is the derivation of estimates not 
involving the singular component of curvature a. 

To illustrate this, we first prove the energy estimates for (3 on H u and for (p, a) on H_ u by 
considering the following set of Bianchi equations: 

V 4 p = div (3 - ■ a + ••• 

V40" = -div *(3 + ■ a + ... 
V 3 /3 = Vp + VV + ... 

However, the curvature component a still appears in the nonlinear terms in these equations. 
In order to deal with this problem, we renormalize p and a. Define 

1 1 

Using the equation 

V 4 X = -a + 
we notice that the first two equations become 

V 4 p = div + V40- = -div *(3 + ... 

At the same time, the equation for /3 can be re-written in terms of of p and a: 

v 3 /3 = Vp + v*o- + ijv^j + ... 

Now we have a set of renormalized Bianchi equations that does not contain a. Using these 
equations, we derive the renormalized energy estimate 

/ /3 2 + / (p,a) 2 < / (3 2 + / (p,a) 2 + / / / V** + V^Wdu'dw', 

JHu Jh„ J H a JH n JO JO Js„, ../ 



in which a does not appear in the error term. 

The same philosophy can be applied for the remaining curvature components (p, a, /?, a). 
As a consequence, we obtain a set of L 2 curvature estimates which do not explicitly couple 
to the singular curvature component a. We say explicitly that there is still a remaining 
possibility that the Ricci coefficients ip appearing in the nonlinear error terms for the energy 
estimates may depend on a. 
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1.2.2. Estimates for the Ricci Coefficients. In order to close the estimates, it is necessary to 
show that all the Ricci coefficients can be estimated without any knowledge of a. The most 
dangerous component is x, which naively would have to be estimated using the transport 
equation 

V 4 X + tixx = -2^X - a. 
We take an alternate route and estimate x from the equation 

V 3 X + ^trxx = V§>?7 + 2ux - ^ X XX + V®V, 

or 

div x = 2 Vtr X - q fa - ZZ) • (X ~ 2 tr *) - 

As we shall see, the loss of information of a is accompanied by the loss of information of 
/3 on H_ u . This presents yet an additional challenge in estimating the Ricci coefficients. 

1.2.3. Higher Order Energy Estimates. Another difficulty arises considering the fact that in 
order to close the energy estimates, we need to prove higher derivative estimates for the 
curvature components. However, the derivatives of some curvature components along e 4 are 
not defined in L 2 initially. We will therefore only use angular covariant derivatives V as com- 
mutators and will prove estimates only for the L 2 norms of the angular covariant derivatives 
of the renormalized curvature components. We will show firstly that this procedure does 
not introduce terms that cannot be estimated in L 2 (in particular a will not appear) and 
secondly that all the energy estimates can be closed only using the estimates of the angular 
derivatives of the renormalized curvature components. 

1.2.4. Existence and Uniqueness. Since we work initial data with very low regularity, the 



a priori estimates that we have derived along the lines of Sections 1.2.1 1.2.2 and 1.2.3 
do not immediately imply the existence and uniqueness of the solutions. Instead, we need 
to approximate the data by a sequence of smooth data and show first that they have a 
common domain of existence and second that they converge in this domain. The fact that 
the solutions to the sequence of smooth data have a common domain of existence follows 
from the a priori estimates that do not involve a and its derivatives as outlined above. Since 
the approximating data are smooth, we can conclude that the approximating solutions are 
also smooth and exist in a common domain. 

Once we have proved the uniform a priori estimates, we proceed to prove that the sequence 
of solutions converges. To this end, we consider the equations for the difference of the 
Ricci coefficients and curvature components. We identify the spacetimes in this sequence by 
the value of their coordinate functions and derive equations for the difference of the Ricci 
coefficients and curvature components. Our a priori estimates heavily relies on the structure 
of the Einstein equations which allows us to eliminate any dependence of the a component. A 
priori, there is no reason to think that this structure is preserved when considering difference 
of these solutions, which is necessary to show that our end result is a strong solution to the 
Einstein equations. It is a remarkable fact that it turns out that in the equations for the 
difference of the Ricci coefficients and curvature components, a indeed does not appear. 
Thus the a priori estimates we have proved are sufficient to control the difference of the 
Ricci coefficients and curvature components from the difference of the initial value. This 
proves that the sequence of solutions converge. Even though we face the standard challenge 
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of loss of derivative in the quasilinear equations, our estimates are still sufficient to show 
that the sequence converges. 

The constructed limiting spacetime is not smooth. In particular, second derivative of the 
metric in the e 4 direction is not even defined. We will, however, show that the Ricci curvature 
tensor is better behaved than a general second derivative of the metric and that the limiting 
spacetime satisfies the vacuum Einstein equations R^ v = in the L? sense. 

The estimates for the difference of the Ricci coefficients and curvature components imply 
that the constructed spacetime is the unique solution to the Einstein equations among the 
class of spacetimes admitting a double null foliation and satisfying strong enough a priori 
bounds. Moreover, we can prove uniqueness of the constructed solution among all limits 
of smooth spacetimes, i.e, any spacetime that arises as a C° limit of smooth solutions to 
the vacuum Einstein equations with initial data converging to the given initial data must 
coincide with the constructed spacetime. 

It can be observed that the above argument for the a priori estimates, as well as that for 
showing existence and uniqueness, does not use the fact that initially, a is a measure whose 
singular support is on So }U ■ In effect, the argument avoids a completely, and can be used 
for data such that a is much rougher. Since the argument used only estimates on the Ricci 
coefficients and the curvature components other than a, it can be used to handle initial data 
satisfying only the assumptions of Theorem [2] 

1.2.5. Regularity and Propagation of Singularity. In the setting of Theorem [TJ i.e., that of 
an impulsive gravitational wave, the theorem gives a precise description of the propagation 
of singularity. The a priori estimates imply that all the curvature components other than 
a are bounded. Here, we are interested in proving two additional statements: firstly, a is a 
measure that indeed has a delta singularity on H u ; secondly, the spacetime is smooth away 
from H„ . 

To show that a is a measure, we approximate the data for a by a sequence of smooth 
data a n such that a n is of size 2 n in a region \u — u s \ < 2~ n . The spacetimes constructed 
for such data are smooth and therefore allow us to use the previously avoided L? estimate 
for the a component of curvature. This estimate imply that in the constructed spacetime, 
in the region \u — u s \ > 2~ n , the L? u norm of a n are uniformly bounded; while in the region 
\u — u s \ < 2~ n , the L? u norm of a n is bounded by 2t. By Cauchy-Schwarz, we have that 
the L\ norms of a n are uniformly bounded. This allows us to show that in the limiting 
spacetime, the curvature component a is a measure. 

To show that the singular part of a is a delta function supported on the null hypersurface 
H„ , we notice that 

a = -V 4 X - tixx ~ 2^X- 

Therefore, it suffices to show that x nas a jump discontinuity across H_ u and smooth every- 
where else. This can be proved using the equation 

V 3 X + ^trxx = V®7] + 2ux - ^t?XX + V®V 

and the fact that on the initial hypersurface H , x has a jump discontinuity across Sq jU and 
smooth everywhere else. 
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In order to show that the spacetime is smooth away from H_ u , we will estimate the higher 
regularity of all the curvature components in that region. We first use the Bianchi equations 

V 3 /3 + trx/3 = Vp + 2uf3 +*Va + 2 X -/3_ + 3(r/p +* rja). 

Integrating this equation and using Gronwall's inequality, we obtain for any u^u s that 

sup^||V73|| L 2 ( ^ M) < (^||V i /3|| L 2 (5o j + ...)exp( / ^||V*(trx 

" i<I i<I J i<I 

The regularity of (3 is thus inherited from the initial data, which is smooth away from u = u s . 
Once we have estimates for /3, we consider the equation 

V 3 a + ^trxa = V®/3 + Aua - 3(xP +* X°) + (C + 4r?)®/3. 

Integrating as before, we see that a also inherits the regularity from the initial data. Higher 
derivatives estimates for a can be derived analogously by differentiating this equation. The 
other components of curvature can be controlled in a similar fashion. Notice that this 
procedure results in a loss derivatives. In particular, in order to control the iV-th derivative 
of the curvature, one needs ~ 2N derivatives initially. 



1.3. Outline of the Paper. Finally, we give the outline of the remainder of the paper. In 
the next Section, we give a careful introduction of the setting, describing the double null 
foliation, the coordinate system, the equations and relevant notations. We will state a more 
precise version of Theorem [2j which we will call Theorem [3j In Section [3j we provide a con- 
struction of the initial data set satisfying the conditions in Theorem [T] and exhibit a sequence 
of smooth data approximating the data with a curvature delta singularity. In Sections [4] and 
[5j we prove Theorem [3] In Section [4], we prove that for smooth initial data satisfying the 
assumptions of Theorem [3j a unique spacetime exists in a region depending only on the 
constants in the assumptions. This in particular implies that the approximating sequence 
of initial data constructed in Section [3] gives rise to spacetimes with a common and uniform 
region of existence (identified by a choice of a double null coordinate system). In Section [5j 
we study the equations for the difference of two spacetimes. This allows us to show conver- 
gence of solutions with converging initial data and conclude the existence part of Theorem 
[3| In Section 5.7 we examine the regularity of the limiting spacetime and show that it is a 
solution to the vacuum Einstein equations. In Section |5.8[ we conclude the uniqueness part 
of Theorem [3] Finally, in Section 5.10 we return to the proof of part (b) of Theorem [Tj 
giving a precise description of the propagation of singularity. 



Acknowledgments: The authors would like to thank Mihalis Dafermos for valuable dis- 
cussions. 



2. Setting and Equations 

Our setting is the characteristic initial value problem with data given on the two charac- 
teristic hypersurfaces H and H intersecting at the sphere S 0:Q . The spacetime will be a 
solution to the Einstein equations constructed in a neighborhood of S^o bounded by the two 
hypersurfaces. 
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Figure 4. The Basic Setup and the Null Frame 



2.1. Double Null Foliation. For a spacetime in a neighborhood of <So,o, we define a double 
null foliation as follows: Let u and u be solutions to the eikonal equation 

gTdpudvU = 0, g^d^ud u u = 0, 

satisfying the initial conditions u = on H and u = on H . Let 

jjp. = -2g^d u u, = -2g^d u u. 

These are null and geodesic vector fields. Define 

2n- 2 = -g(L' : lJ). 

Define 

e% = QL/ , e4 = VtL' 
to be the normalized null pair such that 

g(e 3 , e 4 ) = -2 

and 

l = n 2 iJ, l = n 2 L' 

to be the so-called equivariant vector fields. 

We will denote the level sets of u as H u and the level sets of u and H u . By virtue of the 
eikonal equations, H u and E_ u are null hypersurface. Notice that the sets defined by fixed 
values of (u, u) are 2-spheres. We denote such spheres by S Uj u. They are intersections of the 
hypersurfaces H u and H_ u . The integral flows of L and L respect the foliation S u ^- 

2.2. The Coordinate System. On a spacetime in a neighborhood of S 0j0 , we define a 
coordinate system (u,u, O 1 ^ 2 ) as follows: On the sphere S^o, define a coordinate system 
(O 1 ^ 2 ) for the sphere such that on each coordinate patch the metric 7 is smooth, bounded 
and positive definite. Then we define the coordinates on the initial hypersurfaces by requiring 

—9 A = on H Q , and — 9 A = on H . 
ou ou 

We now define the coordinate system in the spacetime in a neighborhood of So,o by letting 
u and u to be solutions to the eikonal equations: 

grd,nud v u = 0, gTdpvdvU = 0, 

and define 9 1 , 9 2 by 

£ L 9 A = 0, 
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where fii, denote the restriction of the Lie derivative to TS Ut u (See [12] ). Relative to the 
coordinate system, the null pair e 3 and e 4 can be expressed as 

for some b A such that b A = on if , while the metric g takes the form 

g = -2VL 2 (du ®du + du® du) + -f AB (d6 A - b A du) <g> (dO B - b B du). 

2.3. The Equations. We will recast the Einstein equations as a system for Ricci coefficients 
and curvature components associated to a null frame e 3 , e 4 defined above and an orthonormal 
frame ei,e2 tangent to the 2-spheres S Uj u. Using the indices A, B to denote 1,2, we define 
the Ricci coefficients relative to the null fame: 

Xab = g(D A e A ,e B ), X AB = g(D A e 3 , e B ), 

i]A = ~-g(D 3 e A , e 4 ), v A = -^g(D 4 e A , e 3 ) 

w = -^g(D 4 e 3 ,e 4 ), u = -^g(D 3 e 4 ,e 3 ), ^ 

(a = ^g(D A e A , e 3 ) 
where D A = D e . We also introduce the null curvature components, 
a AB = R(e A , e 4 , e B , e 4 ), a AB = R(e A , e 3 , e B , e 3 ), 

/3a = ^R(e A , e 4 , e 3 , e 4 ), §_ A = -R(e A , e 3 , e 3 , e 4 ), 

p = |-R(e 4 , e 3 , e 4 , e 3 ), a = \ *-R( e 4, e 3 , e 4 , e 3 ) 

Here *R denotes the Hodge dual of R. We denote by V the induced covariant derivative 
operator on S u ,u an d by V 3 , V 4 the projections to S UtU of the covariant derivatives D 3 , D4, 
see precise definitions in |22j . 
Observe that, 

u = --V 4 (logfi), u = --V 3 (logfi), ^ 

^ = CA + V A (logfi), ^ = -a + VA(iogn) 

Let ■ (fi^ denote an arbitrary contraction of the tensor product of (jP-^ and (j)^ with 
respect to the metric 7. We also define 

(0 (1) §0 (2) )ab := + 0gVi 2) - <W0 W • (2) ) for one forms $ , , 

A (2) )ab := ^ AB (7 _1 ) C13 0AcA 2 i for symmetric two tensors (/>%, <$ B . 
For totally symmetric tensors, the div and curl operators are defined by the formulas 

(div 4>) Al ... Ar := V B (j) BAl ... Ar , 

(curl ())) Al ... Ar ■.= ft BC V B 4>cA 1 ...A r , 
where j. is the volume form associated to the metric 7. Define also the trace to be 

(tr0) Al ...A r _ a := (7" 1 ) BC 0bca 1 ...a i ._ 1 . 
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We separate the trace and traceless part of x an d X- Let x an d X be the traceless parts 
of x an d X respectively. Then x an d X satisfy the following null structure equations: 



V 4 tr X + ^(tr X ) 2 = -|x| 2 - 2utr X 
V 4 X + trxx = -2^X - a 

V 3 trx+ ^(trx) 2 = -2wtrx- |x| 2 
V 3 X + trx x = ~2uj_X - a 

V 4 trx + ^tr^trx = 2wtrv_ + 2p - X • X + 2div 77 + 2|r Z | 2 ( 5 ) 
1 1 

V 3 tr% + ^trvtrx = 2wtrx + 2p - x • X + 2div r] + 2\r]\ 2 
V 3 X + \^XX = V§7? + 2ux - Jtr^X + 



The other Ricci coefficients satisfy the following null structure equations: 



V 4 r? = -X ■ (V - V) - P 

v s!Z = -x- (v-v) + P 

V 4 w = 2uu + -\r) -77I 2 - ^(77-77) • (77 + 77) - ^|t7 + 77_| 2 + ^p 

3 1 11 

V 3 w = 2ww + -|?7 -77J 2 + 1(77-77) • (77 + 77) - -|77 + 77_| 2 + -p 



The Ricci coefficients also satisfy the following constraint equations 



11 1 
div x = 2 Vtr X ~ ^(V-V)-(X- 2 tr *) ~~ & 

div x = ^Vtrx + ^(77 - 77) • (x - Jtrx) + /? 

2 2 ^ 2 ~ (7) 

curl i] = —curl i] = a + ~x A x 
K = -p + )pc ■ X - ^trxtrx 
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(8) 



with K the Gauss curvature of the surfaces S. The null curvature components satisfy the 
following null Bianchi equations: 

V 3 a + ^trxa = V§/3 + 4ua - 3{xp +* &) + (C + 4rp)§/3, 

V 4 /3 + 2trx/3 = div a - 2u(3 + rja, 

V 3 /3 + trx/3 = Vp + 2u(3 +* Vtx + 2 X ■ (3 + 3(r]p +* r]a), 

V 4 a + hr X a = -div *(3 + l ~x * a - ( ■* (3 - 2 R ■* (3, 
V 3 a + ^trxa = -div + l - X * « - C •* §_ ~ ^ ■* £, 

V 4 p + ^tr X p = div (3 - X -x ■ a + ( ■ (3 + 277 • (3, 

— 3 1 

V 3 P + -tr%p = -div (3 - -x • a + ( ■ §_ - 2r) ■ [3, 

V 4 /3 + trx/3 = -Vp +* Va + 2cof3_ + 2x- f3 - 3(vp -* rp), 
V 3 /3 + 2trxj3_ = -div a - 2uf3_ + 11 ■ a, 

V 4 a + ^trx« = -V<§>/3 + Aua - 3(xp ~* X°) + (C - 4rp)<§/3 

where * denotes the Hodge dual on S UjV ,. 

In the sequel, we will use capital Latin letters A e {1,2} for indices on the spheres S UM 
and Greek letters p G {1, 2, 3, 4} for indices in the whole spacetime. 

In the following it will be useful to apply a schematic notation. We will let <j) denote an 
arbitrary tensorfield, ip a Ricci coefficient and ^ a null curvature component different from 
a. We will simply write ipip (or ip^f, etc.) to denote an arbitrary contraction. The use of this 
notation is reserved for the cases when the precise nature of the contraction is not important 
to the argument. Moreover, when using this schematic notation, we will neglect all constant 
factors. 

2.4. Integration and Norms. Let U be a coordinate patch on So,o and pu be a partition 
of unity in Djj such that pu is supported in Djj. Given a function 0, the integration on S UtM 
is given by the formula: 



/OO rOO 
/ <Ppu ^det 7 d0 1 d0 2 . 
-co J —00 

Let D u#iU by the region < u < tt + , < u < u^. The integration on D Ujjk is given by the 
formula 

p pu pu pco pco 

/ ^'=^2 / ^PuV- det gdd 1 dd 2 dudu 

J Du,u jj Jo JO J -co J -co 

pu pu pco pco 

=2 J2 / ^Pu^V- det lde 1 d9 2 dudu. 

jj Jo Jo J— oo J —oo 

Since there are no canonical volume forms on H u and H_ u , we define integration by 

/' _ rOO rOO 

<t>-=z2 / ^Pu^Vdet^d9 1 de 2 du, 

J H u jj Jo J —oo J —oo 
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and 



H u jj JO J— oo J — oo 



oo roo 



ffiputly/ det -fd9 1 d9 2 du. 



With these definitions of integration, we can define the norms that we will use. Let <f> be 
a tensorfield. For 1 < p < oo, define 



^llz,p(s u ,„) :_ 




<0 

u 


,0>f 


\\<t>\\ P LP(H u ) ■ = 


i H u 


< 0, 


>f , 


101 IWj := 


L 


< 4>, 


>f • 



Define also the L°° norm by 



II0IU~(S^):= snp <0,0>^(0). 

2.5. Precise Statement of the Main Theorem. With the notations introduced in this 
Section, we give a precise version of the statement of Theorem [2] 

Theorem 3. Suppose the initial data set for the characteristic initial value problem is given 
on Hq for < u < and on if for < u < u* such that 

c < | det7 \ Su>0 |, | det7 f 5o , a I < C, 

Ko I + l(4) l 7 \s 0A \)<C, 

i<3 ^ ' 

Oo := ^sup||VV|| i2(5o ,j + ^sup||VV|| i2(5ll , ) < C, 

i<3 - i<3 u 

n « ■= E E supiiv^n^^) + J2 su pii vi *iU 2 (^, ) < c. 

Then for e sufficiently small depending only on C and c, there exists a spacetime (j\A,g) 
endowed with a double null foliation u, u that solves the characteristic initial value problem 
to the vacuum Einstein equations in < u < u*, < u < for u*,^ < e. The metric is 
continuous and takes the form 

g = -2VL 2 (du ®du + du® du) + ^ A B{d9 A - b A du) <g> (dO B - b B du). 

(Ai,g) is a C° limit of smooth solutions to the vacuum Einstein equations and is the unique 
spacetime solving the characteristic initial value problem among all C° limits of smooth 
solutions. Moreover, 

ire 9 ' h> e C ° C ^ L '^ 
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<9 2 <9 2 



dOdu dudu 

In the (u, u, , $ 2 ) coordinates, the Einstein equations are satisfied in L'^L'^L 2 {S). Further- 
more, the higher angular differentiability in the data results in higher angular differentiability. 

We prove this Theorem in two steps. First, we show that if the initial data are smooth 
and the bounds for the initial data in the assumption of Theorem [3] hold, then an e can be 
chosen depending only on the bounds in the assumption of Theorem [3] such that a smooth 
spacetime solving the Einstein equations exists in < u < e and < u < e. This is 
formulated as Theorem |4| and is proved in Section |4j Then, we show that a sequence 
of solutions to the Einstein equations with a converging sequence of smooth initial data, 
satisfying the assumptions of Theorem [3] with uniform constants C and c, converges. The 
limit spacetime satisfies the Einstein equations and has the properties stated in Theorem [3j 
This is formulated as Theorem [6] and is proved in Section |5j We furthermore show that the 
solution is unique and that the regularity in the angular directions persists in Section [5j 

3. The Initial Data 

In this Section, we construct data satisfying the constraint equations such that x has 
a jump discontinuity across a two sphere. We also construct a sequence of smooth data 
satisfying the constraint equations approaching the data with discontinuous X . We derive 
precise bounds for the Ricci coefficients and the curvature components for the initial data in 
this approximating sequence. 

On the initial characteristic hypersurfaces, we fix 

VL = 1 

identically on H$ and H_ Q . This implies 

cu = on H , and u = on if , 

and 

( = r] = —f] on H and H . 
7, x an d X have to satisfy the equations 

£l! = 2 X , (9) 
£l! = 2X, (10) 

^l^x = -\{^x?-\x\ 2 1 (11) 

^Ltrx = -^(trx) 2 -|x|^ (12) 

Here 1^ denotes the restriction of the Lie derivative to TS u ,u. This is a notion intrinsic to 
the null hypersurfaces. 

Prescribe jab, (a, tr% and tr% on the two sphere So,o- Relative to the coordinate system 
(^ 1 ,^ 2 ), we require that on S^q, 

E \(^y(-rAB,(A,trx,trx)\<C. 

i<I+3 

On the initial incoming hypersurface H_ Q) we have the coordinate system (u, l , 2 ) and on 
the initial outgoing hypersurface, we have the coordinate system (u, 1 , 6 2 ) as constructed in 
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Section 2.2 On if , we prescribe smoothly the conformal class of the metric jab satisfying 
y/det jab 1- We require the estimates 

E E k^^-i <- c 

j<J+l i<I+3 

On 5o,o, since jab and 7,4s are in the same conformal class, 

jab = § 2 1ab- 



(10) implies that, 



d 

~q^1ab = 2x AB = 2xab + tr XI AB ■ 

d „9„ d$„ 

— jab = $ ^-Iab + 2^—jab- 
ou ou ou 

{T^falAB = ^logdet 7 = 0, 

—AB = 2 $ dv 1 ^ 

2 <9$ 

depends only on 7 and is in L°°. Thus, (12) can be re-written as 



On the other hand 



Since 



we can identify 



and 



Notice also that 



du 2 8 du du 

This ordinary differential equation can be solved up to u < as long as < e for sufficiently 
small e. The ODE also implies that u < e, $ 7^ 0, and thus the solution is regular. 

On the initial outgoing hypersurface H Q we now prescribe the conformal class of the metric 
jab satisfying yj det jab = 1 such that its u derivative has a jump discontinuity. To this 
end we define smooth matrices (( 7 )i)ab an d ((7)2)^5, such that ((7)1) as is positive definite 
with determinant equals to 1 and both {{j)\)ab and {(7)2) ab satisfy 

E E l(|;)*(|)'(«W«l <- c. 

k<K+li<I+3 — 

Fix Ug < |, where e > is a small parameter depending on the constants C above and will 
be determine later. We prescribe the initial data for 7 such that its u derivative has a jump 
discontinuity across the sphere 5 0jJ4 . Let 

Tab = + (u- u s )((j) 2 )ab1{u>u s }- 

sufficiently small d 

define 



For e sufficiently small depending on C, J AB is positive definite for < u < e. We then 



1 A 
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On 5*0,0, since jab and jab are in the same conformal class, 

lAB = § 2 lAB- 



(|9| implies that 



On the other hand 



Since 



we can identify 



d 

— 7AB = 2XAB = %XAB + tTX^AB- 

ou 

d ^9<S>„ 
—jab = $ -ttIab + 2$— 7 AB . 
ou ou ou 



and 



XAB = l^AB, (13) 

«x = iat (14) 

With this identification, \ has a jump discontinuity at u = u s while tr% is continuous. Notice 
also that 



depends only on 7 and is in L°°. Thus, (11) can be re-written as 



This ordinary differential equation can be solved up to u < u* as long as u^ < e for sufficiently 
small e. (15) also implies that u < e, $ 7^ 0, and thus the metric on the spheres remain 
regular. 

3.1. Approximation Procedure. We now introduce the approximation procedure to con- 
struct a sequence of smooth initial data approaching the data described above. Consider a 
Cq°(WL) function h that is supported in [—1, 1] and is identically 1 in [— |, |]. Let 



2' 2J 



/io(^) — ^o(2s) x > 
x < 



Note that h is smooth. Now let 

n 

j=-oo 

We note that h n is supported in {x > 2~( n+1 )} and /i' n := h n — h n -\ is supported in {2~( n+1 ) < 
x < 2~ n }. Moreover, h n —> l{ x >o} in L p for every p < 00 as n — >• 00 and xh n — > xl^^o} in 
L p for every p < 00 as n -)• 00. At the n-th step, we define 

(7„,Cn,trXn) = (7,C,trx) on 5 , 

and 

(X„,^x) = (X,trx) on H 



20 JONATHAN LUK AND IGOR RODNIANSKI 

as before. Define 

(IJab = ((i)i)ab + (u~ u s )h n (u - u s )((j) 2 )ab on H Q , 

and 

(1u)ab = , ^ (%)ab on H . 
Jdet(j)n 

Notice that for e sufficiently small, we have uniform upper bounds for each component of 7 
and their angular derivatives and a uniform lower bound for det^f. Moreover, it is easy to 
see that for 



7n : = In - Jn-l, 



we have 

"du J y d9 



E !(!;)*( Jj)^»l ^ C2^\ for u s < u < u s + 2~\ k < K + 1, 



?</+3 



Thus 



W fc W^ 1 = °' for - - - s + 2_n ' for * - 1 + 3,k - K + lj 

Moreover, we know that 7„ — > 7 in the sense that 

( f K^) fc (^)'(Tn - 7)l p rfM)" < C2 ( - 1+fc -? )n for A; < K+ 1. 



With the definition of (%)ab, we can solve ( 15 ) for $ n and use (13 ) and ( 14 ) to solve for \ n 
and (trx) n . This then allows us to use the null structure equations and the Bianchi identities 
to solve for all Ricci coefficients and curvature components, as well as their derivatives. In 
the following we will use the notation that a subscript n denotes that metric component, 
Ricci coefficient or curvature component associated to j n . Moreover, we will use the notation 

</>'n = ( t ) n- (Pn-l- 

We now derive the bounds for each of the 0„'s and 0' ra 's in the initial data. By uniqueness 
of the constraint ODEs, all 

4>' n = for u< u s . 
We therefore only derive bounds for u s < u < e. 



From (|15|), we have 

C<*n<C, J2 E K^(^l <^ 
k<K+2i<I+3 



thus 

E E l(^)Utrx)»l<C 



09' 

k<Ki<I+3 



Then using the equation (|15j) again, 

E E l(|:) fe (^)^l<^ 2+ ^ ior u s <u<u s + T 



'du 89' 

k<K+2i<I+3 — 
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and 



>-■• -V$; = o, for u > u s + T n and 2 < k < K + 2, for all i < I + 3. 



du 80 

Then by (fl3| and (Tl4|), we have for x' r 



E E Kl:)^^)'^)^! < ^ fc ", for Ms < « < u s + 2 



'du dO 

k<Ki<I+3 — 



E E K|;) fe (-|nXnM < C2-», for « > u s + 2-; 



and for tr%' r 

r vn 

'du' K d0 



E E ^h k ^ tr ^\ < C2^+*>», for u s < u < u s + 2- 



E E ^(§8^1 ^ C2 ~ n > fOT * ^ 34. + 2 ~" 

fc<ifi</+3 — 

From the equation for ( on Hq together with the Codazzi equation, we have 

1 3 
V 4 C = div x ~ 2 Vtr X - K - 2 tr ^^- 

Notice that since = 1 on H , we can write 



v 4 C = |-C + xC- 

du 



E E i&'&'ci < c. 

fc<X+li<i"+2 



Thus, we can integrate to get the estimate 

+ 1 ~~ 

For the difference (' n = £ n — Cn-i> 

E E i(|)*(^ 4 ci < 

k<K+l i<I+2 — 

Notice that we have used the control on the metric to relate div to the coordinates deriva- 
tives. Since on H , 

C = v = -v, 

we have 

E E <c. 

k<K+li<I+2 — 

and 

E E \(£ l n§j) , K,i)\<c2-"- 

k<K+l i<I+2 — 



By the equation of tr% on H 



1 



V 4 (trx) + -trytrx = 2p - x • X ~ 2div C + Kl 
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and the Gauss equation 



we have 



K = -p + -x ■ X ~ ^trxtrx, 



V 4 (trx) + tr X trx = -IK - 2div C + |Clr 
Since the estimates for the Gauss curvature K can be derived from the estimates for 7 and 
its derivatives alone, by the estimates derived above and the initial estimate for tr% on So,o, 

E E K»sJ * c. 

k<K+li<I+l - 

and 

E E k^'O'Ki s 

k<K+li<I+l - 



By the equation of \ on H 



we have 



and 



1 1 

v 4 x = -2 tr xx - v®c + C®C - 2 tr xx, 

E E s 

fc<_ftT+li</+l — 

E E * 

fc<_ftr+i«</+i - 



By the following equation on H 

V 4 X = -trxx - a, 

we have 

E E i(|-) fc (^) ia »i ^ C2(1+fe) "' for ^ < « < ^ + 2_n > 
E E i(|;)*(^) <a »i ^ C2 ~ n > for * ^ + 2 ~ n > 

fc<K"-li</+3 — 

By the Codazzi equation on H 

1 1 

div x = 2 Vtr X - C • (X - 2 tr x) - ft 

we have 

E E \(§ a n^)%\<Cioru s <u<u s + 2-", 



k<Ki<I+2 

E E \^ k ^e m - C2 ' n for * - ^ + 2 ~ n 



By the Gauss equation on H 



K = -p+ 1 - X .x-\tr X trx, 
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we have 

E E \^m ypnl < C ^u s <u<u s + 2-, 



k<Ki<I+l - 

,d 

'du' y d6 



E E l(|;)*(^)Vnl < C2~ n for u > u s + 2- 



k<Ki<I+l 

By the equation on H 

curl ( = a + -xAx, 

we have 

E E \(^n^n\<Cioru s <u<u s + 2-, 

k<Ki<I+l — 



E E k^) v «i ^ ^2"" for * > «. + 2 ~ n - 

k<K i<I+l 

On ^o^, using the Codazzi equation 



div x = I Vtr X + C • (x ~ Jtrx) + 



2 

we have 

E \(^)%\<ConS 0fl . 

i<I+2 

By the following Bianchi equation on H 

V 4 ^ = -Vp +* Va - tr X P_ + 2x-(3- 3(*(a - (p), 

we have 

E El(|;A|)'4.l^, 

k<K+l i<I — 



By the null structure equation on H 



V 4 w = 3|C| 2 + ^, 



we have 



e tCK^ c - 

k<K+l i=0 — 

y Y\(—) h (—) i J\<C2- 



Together with u = which follows from w = — ^V 4 (logf2), we have the bounds in L°° for 
all the Ricci coefficients and null curvature components on H . 

A similar procedure allows us to solve for the Ricci coefficients and the null curvature 
components on H_ . Recall that for an impulsive gravitational wave, we take the data on H_ 
to be fixed, i.e., ip n = ip and \l/ n = ^ on H_ Q . 
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Finally, we see that on H , we can define V^n an d V^n, for j < J by taking the values 
of V J 3 ip n and V J 3 ^ n restricted to So,o and solve the ODEs as above. Notice that they can be 
solved as long as I is sufficiently large depending on J. 

This allows us to conclude 

Proposition 1. For every I > 2, J, K > 0, there exists characteristic initial data for an 
impulsive gravitational wave such that on H : 

EE E sup ||V|V 4 fc VV|| i2(5o ,j < C ItJ jc, (16) 

j<J k<K i<I+l - 

E E E E su p iiv J 3 v 4 fc v^n i2( ^ o) < c I>J>K , (l?) 

j<J k<K i<I *e{p,cr,/3,a} " 

On Hq, x has a jump discontinuity across So y u s , and 

E E sup ||V J 3 VV|| i2(5o ,j < Cj,j, (18) 
j<j i<i+i - 

EE E supiiviv^iu^^^c^, (19) 

j<j i<i *e{^,p,o-,/3} - 
and for H \ {\u — u s \ > 77} for every r\ > 0, 

EE E supiivMvviu 2( 5o, a )<c w 

j<Jk<Ki<I+l - 

E E E E su p iivjvjv**!^^ < 

j<J k<K i<i te{/3,p,^} - 

and for e sufficiently small, 7 is positive definite on Hq and H^ with bounded angular deriva- 
tives up to the (I + 2)-nd derivative. 

There exists a sequence of of smooth initial data approaching the characteristic initial data 
described above such that |7^)-|7g[) hold for tp n and ^ n , and for H \ {u s < u< u s + 2~ n } ; 

E E E sup ||V|V 4 fc VV n || i2(5o ,j < Cm, (20) 

j<Jk<Ki<I+l - 

E E E E sup II ViVtV^ n | Ua^) < C^ K , (21) 

j<J k<K i<I *e{/3,p,o-,/y- - 

and for H f]{u s <u<u s + 2""} ; 

EE su pII v 3 vj «-IU 2 (^) < Cm? 1 - (22) 

j<j i<i - 

This sequence of data converges in the sense that on H_ Q , the data is fixed; and on Hq, for 
tj}' n =i(> n - i/jn-i and ty' n = ^ n - ^ n -i, 

E E SU P II V 3 V >nl \ms 0& ) < Cj,j2- n , (23) 
j<J i<I+l - 

EE E sup||ViV i <|U 2(5o>Ji )<C / , J 2-« ) (24) 

j<J i<I *e{/3,p,cr,/3} - 
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on H \{u s <u<u s + 2~ n }, 

E E E sup ii^vwjil^ < c IjJiK 2-, (25) 

i<Jfe<Xi</+l - 

E E E E SU P llV J 3V 4 fc V^;|| L2(S0M) < C7,^2-" (26) 

Remark 6. In the context of Theorem [5| a corresponding sequence of ip n and \l/ n can be 
constructed. As opposed to the case of an impulsive gravitational wave, such a sequence will 
satisfy (16)-(19) and (23)-(24), but it will not obey the improved estimates of Proposition^ 



which holds on H \ {u s < u < u s + 2 n }. This additional regularity is only relevant for the 
argument that the impulsive gravitational wave is smooth away from H_ u . 

Proposition [T] in particular constructs the set of initial data satisfying the assumptions 
of Theorem [T] In the following, we will focus on the case I = 2, J = K = 0. This level 
of regularity is sufficient to prove all the estimates and to construct a spacetime metric 
satisfying the Einstein equations. 

4. A Priori Estimates 

In this section, we begin the proof of Theorem [3j We show that for any smooth charac- 
teristic initial data satisfying the estimates in the assumption of Theorem [3j there is an e 
depending only on the constants in the estimates in the assumptions such that the unique 
resulting spacetime admits a double null foliation (u,u) and remains smooth with precise 
estimates on its Ricci coefficients and curvature components in < u < and < u < 
for u,u<e. In particular, no assumptions have to be made on the size of a. 

In the context of Theorem [3j our Theorem in this Section asserts that a region of existence 
can be found independent of n for any smooth approximating sequence of data satisfying 
the bounds in the assumptions of Theorem [3j 

We first define the norms that we will work with. Let 

^ = E E su pii v ^iii 2 (H U )+ E su pII v ^iU 2 (hj 

i<2 We{y8,p,CT >y 9} U *e{p,<r,/3,a} 



O i)P = SUp||VV|| L P (5ll ,„), 



and 



Oi,2= E sup H v VfflU 2 (tf„) + E supIIVV^Hl 2 ^)- 

We write 

2 1 

= E Oi ' 2 + E Oi ' 4 + Oo ' o °- 
The following is the main theorem in this section: 

Theorem 4. Suppose the initial data set for the characteristic initial value problem is smooth 
and satisfies 

c<|det 7 r Su ,J<L7, J2\(h i l\s n , \<C, 

i<3 



2(> 
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C := ^sup||VV|| L 2 (So ti) 

i<2 » 



1Zr 



+ ^sup||VV|| L 2 (Suo) + 

i<2 " 



|V 3 ^|| i2( H ()) + ||V 3 ^|| i2( H )<C, 




IV^I 



L 2 (K ) 



< c. 



*S{p,o-,/3,a} 

Then for e sufficiently small depending only on C and c, there exists a unique spacetime 
(M,g) endowed with a double null foliation (u,u) that solves the Einstein equations in the 
region < u < u*, < u< u^, whenever u*,u^ < e. Moreover, in this region, the following 
norms are bounded above by a constant C depending only on C and c: 

0,6 3 , 2 ,K<C. 

It follows by |36j that for every smooth initial data set, a unique smooth spacetime satis- 
fying the Einstein equations exists in a small neighborhood of Sq^q. In order to show that the 
size of this region only depends on the constants in the assumption of Theorem |4| it suffices 
to establish a priori control in the region < u < e, < u < e. The Theorem then follows 
from a standard last slice argument (See [25]). 

In this section, all estimates will be proved under the following bootstrap assumption: 



0O,oo < A c 



We now outline the steps in proving a priori estimates. In Section 4A 
estimate the metric components under the bootstrap assumption (lA|). In Section 4.2114.5 



(A) 

we will first 



will derive some preliminary estimates and formulae. In Section 4^ we provide Propositions 
which gives L p estimates for general quantities satisfying transport equations. They will be 



we 



used to control the Ricci coefficients and curvature components. In Section 4.3, we establish 



Sobolev Embedding Theorems in our setting. In Section 4.4, we state some formulae for the 
commutators [V4, V] and [V3, V]. They will then be used to obtain higher order estimates for 



general quantities satisfying transport equations. In Section |4.5[ we prove elliptic estimates 
for general quantities obeying Hodge systems. 

After these preliminary estimates, we estimate the Ricci coefficients and the curvature 
components. This proceeds in three steps: 



STEP 1: In Section 
O < C(O ), where C 



4.6 



a 



estimates are proved via the null structure equations 



we prove that 1Z < 00 and 0%^ < 00 together imply that 
is a constant depending only on the initial norm Oq. These 



STEP 2: In Section 4.8, we use elliptic estimates to prove that 1Z < 00 implies 32 < 



C(Oq, TZq, TZ), where C(Oq, TZq, 1Z) is a constant depending both on the initial norm and 1Z. 



STEP 3: In Section 4^, we derive the energy estimates and prove the boundedness of 1Z 
by the initial data. 

These steps provide bounds for the Ricci coefficients up to three angular derivatives and 
the curvature components up to two angular derivatives depending only on the constants 
in the assumptions of Theorem |4j In Section 4.10, we show that this implies that higher 
regularity propagates and that the spacetime remains smooth in a region depending only 
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on the constants in the assumptions of Theorem |4j This allows us to conclude the proof of 
Theorem H] in Section [4. lit 



Finally, in Section 4.12 we return to the sequence of smooth spacetimes arising from the 
sequence of data approximating the data of an impulsive gravitational wave as described 
in Section [3j We show that these spacetimes obey estimates in addition to those given in 



Theorem |4| These bounds will be applied in Section 5.10 to establish the regularity and 



singularity properties of an impulsive gravitational spacetime. 

4.1. Estimates for Metric Components. We first show that we can control Q under the 
bootstrap assumption (lA|): 



Proposition 2. There exists €q = eo(Ao) such that for every e < eo, 

- < n < 2. 

2 ~ ~ 

Proof. Consider the equation 

1 1 X c) 

u> = --v 4 io g fi = -nv.n- 1 = 2^ fi_1 - ( 27 ) 

Fix u. Notice that both u and Q are scalars and therefore the L°° norm is independent of 
the metric. We can integrate equation (27) using the fact that Q^ 1 = 1 on Hq to obtain 



ll^ -1 - l|U~(5 Ul3t ) < C / \\u\\ L °o {Su , jdu < C7A e. 

This implies both the upper and lower bounds for Q for sufficiently small e. □ 
We then show that we can control 7 under the bootstrap assumption dAl): 



Proposition 3. Consider a coordinate patch U on Sq,o and define U Uj q to be a coordinate 
patch on S Uj o given by the one-parameter diffeomorphism generated by L. Define U uu to 
be the image of U u q under the one-parameter diffeomorphism generated by L. Define also 
Du — Uo<«<e o<u<e U U)U . For e small enough depending on initial data and A , there exists 
C and c depending only on initial data such that the following pointwise bounds hold for 7 
in Vjj: 

c < det 7 < C. 

Moreover, in Dy, 

\1ab\Mi- X ) AB \<C. 
Proof. The first variation formula states that 

In coordinates, this means 



From this we derive that 



d 

— 7ab = 2VLxab- 
ou 

d 

— log(det 7) = ntrx- 

}2\ 



Define 7o (m, u, 6\6 2 ) = 7(0, u, 0\ 9 2 

I det 7 - det (70) I <C \tr X \du! < CA e. 

Jo 
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This implies that the det7 is bounded above and below. Let A be the larger eigenvalue of 
7. Clearly, 

A < C sup 7, (28) 

A,B=l,2 



and 

^2 \Xab\ 2 < CA\\x\\l° 

A,B=1,2 

Then 



nil 

\1ab ~ (1o)ab\ < C I \xab\(1u < CAA e. 
Jo 



Using the upper bound (28), we thus obtain the upper bound for |tab|- The upper bound 
for |(7 _1 ) AB | follows from the upper bound for |tab| and the lower bound for det7. □ 

A consequence of the previous Proposition is an estimate on the surface area of each two 
sphere S u& . 

Proposition 4. 

sup \Area(S Ut u) - Area(S ufl )\ < CA e. 



Proof. This follows from the fact that -y/det 7 is pointwise only slightly perturbed if e is 
chosen to be appropriately small. □ 

With the estimate on the volume form, we can now show that the LP norms defined with 
respect to the metric and the LP norms defined with respect to the coordinate system are 
equivalent. 

Proposition 5. Given a covariant tensor $a x ...a t °n <3u,u, we have 



II ^...AlV^fde'de 2 

i=l Ai=l,2 J "* 



We can also control b under the bootstrap assumption, thus controlling the full spacetime 
metric: 

Proposition 6. In the coordinate system (u,u,6 l ,8 2 ), 

\b A \ < CA e. 

Proof. b A satisfies the equation 



du 

This can be derived from 



f)h A 

---4n 2 C- 



ILM 9bA 8 



du 06 A 

Now, integrating and using Proposition [5] gives the result. □ 
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4.2. Estimates for Transport Equations. The estimates for the Ricci coefficients and 
the null curvature components are derived from the null structure equations and the null 
Bianchi equations respectively. In order to use the equations, we need a way to obtain 
estimates from the null transport type equations. Such estimates require the boundedness 
of trx and trx, which is consistent with our bootstrap assumption (|Aj). Below, we state two 
Propositions which provide L p estimates for general quantities satisfying transport equations 
either in the e% or e 4 direction. 

Proposition 7. There exists eo = eo(Ao) such that for all e < eo and for every 2 < p < oo, 
we have 



IM|i*(S Uia ) < C (\\4>\\lp(s uu ,) + / \\^4(p\\LP(S uu ,,)du') 



LP 



Proof. The following identity holds for any scalar /: 



d 

du 



f = f (£ + Otrx/) = / n (e 4 (/) + trx/) • 



Similarly, we have 



/ /= / 0(es(/)+trx/). 

" SlL.1l, J SlL.1t 



Hence, taking / = we have 



LP(S u>!k ) — IIVllii»(S ,) 



+ / / p|0r 2 fi (< 0, V 4 > 7 +-trx|0|^ rfn" 
P \(Pr 2 n (< 0, V 3 > 7 +^trx>|f) du" 



(29) 



The Proposition is proved using Cauchy-Schwarz on the sphere and the L°° bounds for 
f2 and trx (t r x) which are provided by Proposition [2] and the bootstrap assumption (|A 
respectively. 

The above estimates also hold for p = oo: 
Proposition 8. There exists eo = eo(Ao) such that for all e < eo, we have 



L-(5„, a ) < C ( 1 10| | l~ (5^,) + y ^ ||V 4 0|Uoo ( c ?u ^ / , ) cZw ,/ 

L-( Su , M ) < c M|0IU~(s a ,,j + y ||v 3 0|U- (Sii „ jdw"^ . 

Proof. This follows simply from integrating along the integral curves of L and L, and the 
estimate on Q in Proposition |2j □ 
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4.3. Sobolev Embedding. Under the bootstrap assumption (|Aj), the Sobolev Embedding 
hold on a 2-sphere S UtU . 

Proposition 9. There exists €q = eo(Ao) such that as long as e < €q, we have 



i 



L*(Su,u) < °Z^H N ' ; Nl~!.s' ; , LL ,. 
i=0 



Proof. We first prove this for scalars. Given our coordinate system, the desired Sobolev 
Embedding follows from standard Sobolev Embedding Theorems and the lower and upper 
bound of the volume form. Thus, the proposition holds for scalars. Now, for <fi being a 

tensor, let / = \l\4>\^ + 5 2 . Then 



L*{S UA ) < C I ll/IU 2 (S u , 2 ) + ||— F^=^=^IU 2 (Sn, a ) I < C (ll/IU 2 ^) + l|V0||L2 (5tiM) ) 



The Proposition can be achieved by sending 8 — > 0. □ 

We can also prove the Sobolev Embedding Theorem for the L°° norm: 
Proposition 10. There exists to = eo(A ) such that as long as e < eq, we have 

U\\l^(s ua ) < C (\\<f>\\L*(s v , % ) + ||V0|| L 4 (Sufi) ) . 

As a consequence, 

2 
i=0 

Proof. The first statement follows from coordinate considerations as in Proposition [9j The 
second statement follows from applying the first and Proposition [9j □ 

4.4. Commutation Formulae. We have the following formula from [22] : 
Proposition 11. The commutator [V4, V] acting on an (0,r) S-tensor is given by 
[V 4 , Vb]^..^ =[D4,D B ]<l> AlmAr + (V B logfi)V 4 0A 1 ...^ r - (7~ 1 ) CD XB D Vc<f>A 1 ...A r 

r r 

- J2^ CD ^DV Ai ((> Al _„ Ai a„ Ar + J2^~ 1 ) CD XA i BV D (l> Al .„ Aia ,. Ar . 

i=l i=l 

By induction, we get the following schematic formula for repeated commutations: 
Proposition 12. Suppose V 4 = F . Let V 4 V 4 = Fj. Then 

il+«2+«3+i4=*-l 

where by V ll (f] + ?7) i2 we mean i/ie stzm of all terms which is a product of 22 factors, each 
factor being V J (r/ + rf) for some j and that the sum of all j 's is i\, i.e., V ll (n + rf) 12 = 



Jl + ...+j'io=»l 
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V jl (r] + 77). ..V^ 2 (rj + n). Similarly, suppose V30 = Go- Let VaV*^ = Gj. Then 

Gi ~ ^ V il (^ + !/) 42 V i3 Go+ J] V il (?; + 2) <a V <8 xV <4 
+ + # V* 3 /3V* 4 0. 

«l+«2+«3+*4=i— 1 

Proof. The proof is by induction. We will prove it for the first statement. The proof of 
second one is analogous. The formula obviously holds for i — 0. Assume that the statement 
is true for i < i . 

F t0 =[V 4 ,V]V i °-V + VF i0 _i 

~y\>0 + (77 + 5) V 4 V i0 - V + /3V io ~V + x(r? + rj)V i0 - l <j) 

«l+«2+«3=io 

+ V il (?; + 2)* a V <3 xV i4 

«l+«2+*3+M=«0 

+ V^ + # 2 V 43 /3V i4 0. 

11+212+13+14 =*o— 1 

First notice that the first, third and fourth terms are acceptable. Then plug in the formula 
for Fj _i = V 4 V io_1 an d get the result. □ 

The following further simplified version is useful for our estimates in the next section: 
Proposition 13. Suppose V 4 = F . Let V 4 V*0 = Fj. Then 

Fi~ V il # 2 V i3 F + V^^V^xV* 4 ^ 

H+i2+i3=« ii+*2+*3+H=i 

Similarly, suppose V30 = Go- V3VV = Gj. T/ien 

Gi~ ^ V il # 2 V i3 G + V^^V^xV* 4 ^. 

il+«2+i3=« U+i2+*3+«4=i 

Proof. We replace /3 and /3 using the Codazzi equations, which schematically looks like 

= Vx + ^X, 
13 = Vx + Vx- 

□ 

4.5. General Elliptic Estimates for Hodge Systems. We recall the definition of the 
divergence and curl of a symmetric covariant tensor of an arbitrary rank: 

(div 0)Ai...Ar = V B (j)BA 1 ...A r , 
(curl 4>)A 1 ...A r = i BC y B<PcA x ...A r , 

where £ is the volume form associated to the metric 7. Recall also that the trace is defined 
to be 

The following elliptic estimate is standard (See for example p2] or [12]): 
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Proposition 14. Let (f) be a totally symmetric r + 1 covariant tensorfield on a 2-sphere 
(S 2 ,7) satisfying 

div (p = f, curl <f) = g, trcj) = h. 

Suppose also that 

J2\\^ iK \\L^S) < OO. 
i<l 

Then for i < 3, 

i-l 

llVVIUs(S) < C(J2 l|V fc ^|U2 (s) )(^(||VV|| L2(s) + ||V^|U 2(5) + ||V^|| L2(S) + 

k<l j=0 

For the special case that <fi a symmetric traceless 2-tensor, we only need to know its 
divergence: 

Proposition 15. Suppose <fi is a symmetric traceless 2-tensor satisfying 

div 4> = f. 

Suppose moreover that 

^llV^H^s) < oo. 

i<l 

Then, for i < 3, 

i-l 

llvviUs(s) < c(£ \\v k K\\ LHS) )(Y^(\\v'f\\ LHS ) + \\<P\\ms)))- 

k<l j=0 



Proof. In view of Proposition 14 , this Proposition follows from 

curl =* /. 

This is a direct computation using that fact that <fi is both symmetric and traceless. □ 

4.6. L P (S) Estimates for the Ricci Coefficients. We continue to work under the boot- 
strap assumptions ([A]). In this Subsection, we pursue step 1 of the proof as outlined in the 
beginning of the Section. We show if that the curvature norm 1Z and the third order Ricci 
coefficient norm 0^,2 are bounded, then so is the Ricci coefficient norm O. In particular, our 
bootstrap assumption (|A]) and all the estimates in the last section are verified as long as 1Z 
and 03,2 are controlled. 

Proposition 16. Suppose 

1Z < OO, 03 ; 2 < OO. 

Then there exists e = e o (0o, ^3,2) such that for e < e , 

O < C{O ), 

where Oq is the initial norm for the Ricci coefficients defined in Theorem^ 

Proof. To prove the Proposition, we need to estimate 0^4 and O^- First, we estimate 

i<l i<2 

the Oi t 4 norms. Using the null structure equations for all Ricci coefficients except x, we have 
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and 

Y V 4^= Y *+^, 
ipe{ti x ,ri,ui} ve{p,p,(r,P} 

where the ip on the right hand sides denote any Ricci coefficients. We will also use the null 
structure equations commuted with angular derivatives: 

V^= Y Y V n ^ i2 V i3 ^+ V V^V^V* 4 ^, 

t/>e{trv v T) ujI *e{p,tr,/3,a} U+«2+i3=* n+«2+i3+«4=« 

- (30) 

jp£{tlx,ri,ui} *G{/3,p,o-,/3} ii+«2+«3=« ?i+i2+«3+*4=« 



where we have applied the notation defined in Proposition 12 By Sobolev Embedding in 
Proposition [9] and the bootstrap assumption (LAI), we have 



Y HV 4l ^ 2 V i3 vl/|| L4( ^ M) < A ^||V^|U 2( 5 Ui£ ). 

U+i2+«3<l 1=0 

1 

Y \\V^V^V^\\ LHSu ^ < (A + Al)YO iA . (31) 

*l+«2+i3+*4<l 1=0 

Hence, by Proposition [7] applied to plj ), we obtain 

1 

Y ®iA tr X, V, V, u, w, tax x] 



i=0 

1 

<O + Ce^A TZ + C7(A )e £ C 4 



(32) 



'i,4- 
i=0 

In order to close this estimate, we need to estimate x- We like to avoid the use of the 
equation V4X = —a + ipip since we do not have an estimate for singular curvature component 
a. Instead, we use the equation 

Commuting with angular derivatives, we get 

V 3 V i x= Y V il ^ 2 V i3+1 ?7+ Y V^V^V* 4 ^- (33) 



Applying Proposition 7 to (33) and using the estimates in (31), we have 



EllV**!!*^ ^ C ^ \\^ 2 v\\L H s uf jdu' + C(A )eYO tA . 

i=0 ^° i=0 

By Cauchy-Schwarz and the Sobolev Embedding in Proposition [9j 

/ \\V 2 v\\lhs u , u )du' < Ce0 2 , 2 + Ce 1 z6 3 , 2 . 
Jo 

Therefore, 

1 1 
Y H V ^IU 4 (^) < C ^°2,2 + Ce5(9 3i2 + C(A o )e^0 M . (34) 

i=0 i=0 
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Putting (32) and (34) together, we get 



°^ + C^AoTl + Ce0 2 ,2 + Ce50 3 ,2 + C(A )e ^ O iA . 



(35) 



i=0 



i=0 



By choosing e sufficiently small depending on Aq, Oq, ^ and 03,2, we have 



Y,O iA <C(0 Q ) + CeOw, 

i=0 

where C(Oq) does not depend on A , once e is chosen to be small enough. 

We now move to estimate the 0j2 norms. Applying Proposition [7| to (30), we obtain 



^ i tT X, V, V, w, to, trx, x] 

i=0 

2 1 
<O + A TZ + eA ®i,2 + <Y1 

i=0 i=0 

For x, we apply Proposition [7] to (33) to get 

2 21 
E W^xWlHS^) < C(O ) + Ce^6 3 , 2 + CeA ^ (9 i>2 + e(£ 



(36) 



(37) 



i=Q 



i=0 



i=0 



(36) and (37) together give 

2 21 
J] O h2 <O + Ce^A TZ + eA £ 4 , 2 + e(^ £\ 4 ) 2 - 

i=0 i=0 i=Q 

Choosing e to be appropriately small depending on A , O , 1Z and 03,2, we can absorb 
2 

eA Oj,2 to the left hand side to obtain 



i=Q 



(38) 



i=0 



i=0 



For e sufficiently small, (35) and (38) imply that 

1 2 



i=0 



8=0 



Finally, we use Sobolev embedding in Proposition 10 to get the necessary bounds for Oo )OC 
1 

from the control for Oj,A- Moreover, Ao can be chosen so that Ao ^> C{Oq) and the 

i=0 

bootstrap assumption |a]) is verified as long as 1Z and O^^ are controlled. Furthermore, the 
since the choice of A depend only on O , the e can be set to depend only on O , 71 and 



3,2- 



□ 
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4.7. L 2 (S) Estimates for Curvature. In this Subsection, we prove that f3,p,a,f3 and 
their first derivatives can be controlled in L 2 (S). The argument relies on Proposition]?] and 
the Bianchi equations pi). The results of this Section will give the bounds for the Gauss 



curvature and its first derivatives appearing in Proposition 14 
Proposition 17. Suppose 

1Z < oo, 3j2 < oo. 
Then there exists eo = eo(Oo, ^3,2) such that for e < e$, 

E E \\v i nLHs u , uJ <c(n ). 

Proof. By the Bianchi equations, we have 

v 3 G0,p,<7,£)= E w + ^ E * 

Commuting the above with V, we have 

V 3 V(Ap,<7,£)= E y2vI/+ E (*V^ + ^V*). 

Estimating directly, we get 

El|V 3 VX/3,p,a,£)|| x2(s j 

<E E iiv^iu^j+csupEiiv^iu^e E 11^*11^) 

i<2 fe{p,o-,/3,a} " ii<l «2<2 *e{p,o-,/3,a} 

+ Ce sup E IIV^IU*^) E E 1 1 V " *l I^CAmJ- 

" ii<l i 2 <2 *e{p, CT ,/3,a} 



Therefore, by the definition of 7?. and Proposition 16, we have 
1 

E 1 1 V 3 V*(/3, p, a, P) \\ LHK j < C(O )TZ + C(O )e E 1 1 V l (/3, p, a, £) \ \ L , {s ^y 

i=0 i<l 

Therefore, by Proposition [7j we have 

E 1 1 V*(/3, p, a, £) I \ LHSu ^ < C(K ) + e^C(O )7l + C(O )el E 1 1 V l (A P, a, £) | \ L * {s ^j. 

i<l i<l 

The conclusion follows from choosing e sufficiently small and absorbing the last term to the 
left hand side. □ 

We now derive L 2 (S) estimate for K and its derivative, which will in particular allow us 



to apply Proposition 14 



Proposition 18. Suppose 

1Z < oo, (9 3>2 < oo. 
Then there exists e = eo(Oo, 71, 0^,2) such that for e < e , 

Ellv^lMs) <c(o ,n ) 



i<l 



36 



JONATHAN LUK AND IGOR RODNIANSKI 



Proof. We use the equation 



K = -p + ■ X - j-trx ■ tT X, 



and bound term by term using Propositions 16 and 17 



□ 



4.8. Elliptic Estimates for Ricci Coefficients. This Subsection contains step 2 of the 
proof of the a priori estimates outlined in the beginning of the Section. We bound the third 
derivatives of the Ricci coefficients. They cannot be estimated by transport equations alone 
as before because that would require the control for the three derivatives of curvature, which 
we do not have. Instead, we need to use elliptic estimates, following ideas in [T3], [22], [12], 
[21]. We hope to find carefully chosen combination of G := Vip + \1/ such that satisfies 
transport equations 

V 3 6 = ipVipH + iH>*/> + i/>*h, (39) 

or 

v 4 e = ipVip H +ii)ipip+ip^H, (40) 

where ip H G {tr*, x,rj, % u, trx), ipH G {trx,V,V,^.,^X,x}, G {/3,p,a,/3} and ^ K e 
{p,a, (3, a}. These transport equations will then be coupled to the estimates from Hodge 
systems of the type 

Vtp = e + ^ + ipip 

to recover the control for the third derivatives of the Ricci coefficients. In order to apply 



this strategy, it is desirable for (39) and (40) to have two special structures. Firstly, V\l/ 



should not appear on the right hand side of the equation so that can be estimated by 
these transport equations without a loss of derivative. Secondly, we hope that the terms 
ifjVipH and ip^H do not enter the V30 equation (respectively, ipVipH_ and ip^H_ do not enter 
the V40 equation). If this is satisfied for equations p9|) and (40), then the integrated error 



terms can be controlled by the flux term that are bounded on the corresponding hypersurface. 
Indeed, the third derivatives of most the Ricci coefficients can be obtained in this way by 



constructing approximate variables which satisfy either (39) or (40). 



An extra difficulty arises when we estimate V w, More precisely, the variable that 
allows us to retrieve estimates for derivatives of u satisfies an equation of the form 

V 3 e = if)VipH+iJnW + i/)VH + 'ipP, (41) 

which has an undesirable term ipfi. Since we prove estimates without any information on the 
singular curvature component a, we automatically lose control of V 2 /3 in L 2 (H_). Therefore, 
when obtaining the estimates from this transport equation using Proposition [7j the error 
term arising from ip/3 in (41) cannot be bounded. In order to overcome this problem, we 
prove a weaker bound for V^w. More precisely, we will only control it in L 2 (H) but will not 
obtain any bounds of its L 2 (H_) norm. We will moreover show that this is sufficient to close 
the energy estimates. 



We first bound the variables that satisfy equations of the form (39) and (40) 



Proposition 19. Let satisfy the equation 

V 3 = 1pV^) H + Iplplp + tp^H, 



or 



V 4 = TpVlpH + IpTplp + Ip^H, 
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where ip H E {trx, X, V, V, trx), 4>h e {trx, rj, 5, u, trx, x), G {0,p,a,0} and ^ K G 
{p, a, (3, a}. Suppose 

iz < 00, d) 3i 2 < 00. 

Then there exists e = e (O , 1Z, 3t2 ) such that for e < e . 

2 

2=0 — 

Proof. Notice that the assumption allows us to choose eo such that the conclusion of Propo- 
sition 16 holds. For satisfying 

V 3 = ijNipH + irW> + ^K, 
we commute with angular derivatives to obtain 

»X+»2-H3-H4=» 

Take % < 2. We estimate term by term. For the first two terms, except for the factor with 
V 3 i/jh_, can be estimated using the O norm, which is bounded in view of Proposition 16 The 
term ipV^ipn can be controlled by the norm. We have 

£ HV^V^V^+Vh + V il ^ 2+1 V l ¥V i4 ^|| L i L2( 5 tl „ ) < C(O )(l + ^63,2). 

«l+«2+i3+*4<2 



The curvature term can be estimated using Proposition 16 and the definition of the curvature 
norm, since the component f3 does not appear: 

llV^'V^V^I^s^) <C(O )(l + e^TZ). 

«l+*2+«3+*4<2 

We then estimate using Proposition [7j 

2 

sup liv'eiu*^) < c(o ) + c(o )ekn + 6 3j2 ). 

1 = 

Similarly, for satisfying 

V 4 = IpVlpH + lH>$ + ^H, 

we commute with angular derivatives to obtain 

V 3 V*0 = £ (V^^V^V^+Vtf + V <1 V' <a+1 V <s VV i4 V + V h ip h V i3 ^V h ^ H ), 

We then estimate as above and use Proposition [7] to get, 
2 

£ sup ||v*e|| L 2 ( s ttiw) < c(o ) + c(o )ekn + o 3)2 ). 

■ n n 

i=0 — 

The Proposition follows from choosing e sufficiently small depending on IZ and O^^- D 

We now write down the explicit quantities following [21] and show that the estimates 
for some of the Ricci coefficients can be recovered from via elliptic equations. 
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Proposition 20. Assume 

1Z < OO, &3 t 2 < OO. 

Let G denote any of the following quantities: 

Vtrx, Vir%, p := —div r\ — p, p := —div 77 — p. 



Then 6 satisfies (39) or (40). Moreover, there exists eo = e (Oo,Tl,03 : 2) such that for 
e < eo? we have 

sup||V 3 (irx, trx)\\ms^) < C{O )eS, 



u,u 

73/ 



sup ||V J (x,^7,ZZ)IU 2 (^) < C{Oo)(e* + 

u 

sup ||V 3 (x,^,!7)IU 2 (j/ u ) < C(O )(e*+K), 

u 

where, as before, ip H G {trx,X,V,V,u, trx), fee {trx,V,V,u, trx,x}- 
Proof. Consider the following equations: 

V 3 trv_ = ipHipH, 

V 4 trx = ipHipH- 

These equations do not involve a curvature term. After commuting with an angular deriva- 
tive, we have 

V 3 Vtrx = ipVipH + ipipip, 

V 4 Vtr% = ipVipH + 4>4>4>- 

Thus Vtrx and Vtr% satisfy the requirement to be one of the G quantities. Thus, by 
definition, 

l|v 3 trxlU 2 (5 u , a )<c||v 2 e|U 2(5uia) . 

Moreover, 

||V 3 tr X || i2(ifn) <Ce5||V 2 e|U 2(5ttia) . (42) 

and 

l|V 3 trx|U 2( ^j<Ce^||V 2 eiU 2( ^ a) . (43) 

From the estimates for V 3 tr% and V 3 tr\;, we can obtain the estimates for V 3 x and V 3 %. 
This is because of the Codazzi equations: 

div x = ^Vtrx - (3 + ipip, 
div x = 2 Vtr X + §_ + 

Since x an d X are symmetric traceless 2-tensors and we have the estimates for the Gauss 
curvature in Proposition 18, we can apply the elliptic estimates for symmetric traceless 
2-tensors in Proposition [15] 

t3 j ll ,ll t — t 1 ?. r\ \ l X l l t — ri-\ i t — Tin i l l \ 



|V 3 x||l 2( ^) <C(||V 3 tr X || L2(Hii) + ||V 2 /3|U 2{ ^)+ W^VHWl^ 



U+i 2 <2 
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Since, by Proposition [16 



IIv^v j ¥IIl W 

n+ T~ 2 , (44) 

we have 

I |V 3 x| \l*(h u ) < Ce l 2 1 1 V 2 9| U 2(5uia) + CTZ + C(O )ei . (45) 
Similarly, using the Codazzi equation for x, we get 

||V 3 x|| L2(i y < Ce5||V 2 0|| L 2 (Su ^) + CTZ + C(O )e^ . (46) 

Define 

p = — div r] — p, 
p = —div T] — p, 

where as before p = p — \x ' X- Notice that div rj and div rj satisfy the following equations: 

V4div rj = —div (3 + ipVip + ipipip + ip^H, 

V 3 div rj = div (3 + ipVip + ipipip + ip'Q h_, 
Recall the Bianchi equation: 

V 4 p + ^tr X p = div - X -x ■ a + ( ■ P + 2 R ■ (3, 

3 1 
V 3 P + 2 tr XP = -div p - -x ■ a + ( ■ p - 2rj ■ p, 

In the V4 equation p, the term a appears. This is undesirable because we do not have good 
estimates for a. Instead we renormalize the equations using 

1 . A 
P'=P- ^x-x- 

Now, we write down the equations for p instead of a and p. Since V4X = —a + ipip, the V4 
equations now do not contain a. Schematically, we have 

V 4 p = div p + ipVip + tp^H + ip 3 , 

V 3 p = -div p + ipVip + ip^H + V' 3 , 

where all the curvature components are now acceptable. Thus p and p satisfy the following 
equations: 

V4P = Ip^H + Ip^lpH + 4>ll>1p, 

V 3 p = + ^V^Ph + 

Thus p and p also satisfy the schematic equation for 0. We now show that we can get the 
estimates V 3 // and V 3 ?? from that of V 2 p and V 3 p. To do this, we consider the following 
div-curl systems 

div rj = — p — p + ipip, 

curl rj = a + ipip, 
div 77 = — p — p + ipip, 
curl rj = —a + if}tp. 
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Notice that p and a can be estimated by 1Z on both H u and H_ u . Then, using (44) and 

(47) 



Proposition 14, we have 



and 



\^(v,v)\\lhh u) < Ce5||V 2 ei| i2(5 „, H) + cn + C(O )e^ 
|V 3 (7 7)!? )|| L2( ^ ) <Cel||V 2 e|U 2(5uj3i) + C'^ + C(C»o)ei 



(48) 

The Proposition thus follows from (|42), (|43j), (|45j), (|46j), (|47j), (|48j) and Proposition [l9j □ 

The above construction of provides the desired estimates for the third derivatives of all 
the Ricci coefficients with the exception of u and u. Therefore, to estimate the O^^ norm, 
it remains to prove estimates for V 3 u and V 3 cu. To this end we follow the strategy in 
Define two auxiliary quantities and ur to be solutions to 

V 3 w t ' 



V 4 w f 



2 ' 



with zero boundary conditions along Ho and H_ respectively. For the remainder of the 
Section, ip will also be used to denote and u}; ipH will also be used to denote a/; and ipn 
will also be used to denote uj} . It is easy to see that 

Proposition 21. Suppose 

K<oo. 

Then there exists e = e (Oo,7V) such that for every e < e , 

^\\V(u\J)\\ws^)<C. 



i<2 



Proof. Since 



V 4 w t 



2 (7 ' 



:cr. 



with zero data. It follows from Proposition [7] that 

^||V^a;t,a;t)|| i2(Suia) <c(Oo)6^. 

i<2 



□ 



The Proposition follows from choosing e appropriately. 

Now are now ready to use elliptic estimates to control V 3 u and V 3 o;. As before, our 
strategy is to use a combination 9 = Vip + \P that can be estimated using a transport 
equation. However, in this case does not satisfy an equation of the type in Proposition 

us 

Proposition 22. Suppose 

1Z < OO, 0$ t 2 < oo. 

Then there exists eo = €q(Oo,TZ, 03,2) such that whenever e < eq, we have 

\\V 3 (u,^)\\ LHHu) < C(O )(e^ +n), 
\\V 3 (u,(J)\\ lH hj < CiOo)^ +11), 
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Proof. Let 



We have 



and the Bianchi equation 

Therefore 

Similar, using 



and the Bianchi equation 
we have 



re := Vw +* VuA - 

V 3 Vu; = ^Vp + VW> + V* + ipipip, 
V 3 * Vw f = ^*V(x + ^ + ^Vu f , 



Vp +* Vcr = V 3 /3 + V*- 
V 3 re = + V'W + iiipip- 

v 4 Vw = ^Vp + vw> + + 

V 4 *Vw t = ^*V(x + ^ + ^Vu f , 
-Vp +* Vo- = V 4 /3 + v*, 



(49) 



(50) 



V 3 re = ip'fy + ^V^ + 

The equations for re and re do not have the favorable structure of the equations of the other 
6 variables. In particular, for the V 3 re equation, (3 and Vu appear as source terms. The 
loss of the L 2 energy estimate for a is automatically accompanied by the loss of estimate for 
(3 in L 2 (H). Therefore, V 2 /3 cannot be controlled in L 2 (I£ U ). This means that we cannot 
estimate V 2 re in L 2 (H_ U ). Nevertheless, as we now show, we can estimate V 2 re in L 2 (H U ). 
This estimate will be sufficient in the sequel. 
Commuting (49) with angular derivatives to get 

V 3 V 2 re= ^ n V i2 ^V i3 ^+ Y ^V^V^. 

il+«2+i3<2 U+»2+23<3 



Notice that we now allow ip to represent the terms and ur . Since by Proposition 21 
and V 1 uj} satisfy the same estimates as any other ip for i < 2, those terms can be estimated 
in the same way as in Proposition [19} The new terms that have not been estimated in 
Proposition 

H+l2+i3<2 



are ipV 3 uj} , ip n V t2 ipV t3 (3 and ipV 3 (x, u, w*). We estimate them by 



^VV||l!l2 ( s) < C(O )e3||VVIU^ 



(S), 



and 



and 



(5Vi,;u) ' 



i<2 



||^V 3 (x,^)|| L i L2(s) < C(0 )e>\\V 3 (x,oj)\\ LlL2 



(s)- 
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Then, using Proposition [7] and the estimates in the proof of Proposition 19, we have 



llv^lU^ 



<C(O )(l + ^03,2 + e*1l + ||VV|| L 2 L2(5um) 

i<2 

This holds for any fixed u, u. Integrating this in L 2 , we get 

\\v 2 K\\mH u) 

<C(O )e^ (1 + e5 (a 3>2 + ft + || VV| 

+ ^l|V l /3||L i L^(5) + l|V 3 (x,w,u; t )|| LiL 2 L2(5) ). 

Using Cauchy-Schwarz in the w variable, we have 

j<2 

and 

Therefore, 

\\V 2 k\\ L 2 {Hu) < 0(0^(1 + e*d 3 , 2 + e^ + e 5||VV||L r L^ 2 (^)+^||VV|| i ^ L 2 i2(5iiiti) ). 

(51) 

Similarly, we can use (50) to conclude that 

1 1 V 2 «:| \ L 2 {K j < C(0 o )ei (1 + (5 3)2 + + \ | V V | |^ig^(s„ >3i ) +^|| WlU^a^)). 

(52) 

We show that we can recover V 3 w, V 3 w, V 3 wt and V 3 wJ" by considering the div-curl system 
for Vw and VuA We have 

div Vw = div k + T^iv ^' 
curl Vw = 0, 
div Vw^ = curl k + -curl (3, 

curl Vw^ = 0, 
div Vw = —div k — 0, 

curl Vw = 0, 
div Vw = —curl k — -curl (3, 
curl Vw = 0. 



Thus, by Proposition [14] and 18 we have 

||V 3 (w,wt)|| i2( ^ ) <C(||V 2 K || LW + ^||V l /3|| i2( ^ ) ), 



i<2 
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and 

liv 3 ^)!!^) < c{\\v\\\ L , {K j + £ ||v^|U 2( hj). 

i<2 



Using (51 ) and (52), we can thus choose e sufficiently small to conclude the Proposition. □ 



Putting these all together gives 
Proposition 23. Assume 

n<oo. 

Then there exists €q = (Oo,7Z) such that 

o 3 , 2 <c(o )^+n). 

Proof. Impose the bootstrap assumption 

03,2 < Al (53) 
Under this assumption, we can choose e depending on C , TZ and Ai such that the conclusions 



in Propositions 19 20 and 22 hold. Thus, 

d 3 , 2 <c(o )^+n). 



This improves the bootstrap assumption (53) for Ai sufficiently large. Moreover, Ai can be 



chosen to depend only on Oq and 1Z. Thus, eo depends only on Oq and 1Z. □ 

4.9. Energy Estimates for Curvature. In this Subsection, we carry out step 3 in the 
proof of the a priori estimates as outlined in the beginning of the Section. We prove en- 
ergy estimates for the curvature components and their first two derivatives, thus providing 
bounds for 7Z and closing all the estimates. These energy estimates will be derived after 
renormalization directly from the Bianchi equations. This is analogous to [TS] , where a 
similar procedure without curvature renormalization has been carried out to derive energy 
estimates for all components of curvature. We remind the reader of our challenge to avoid 
any information on the non-L 2 -integrable curvature component a in our context. 

We recall the following two Propositions from |28j, which can be proved by a direct com- 
putation. 

Proposition 24. Suppose 0i and 02 are r tensorfields, then 

0lV 4 02 + / 02V 4 01 = / 0102- / 0102+ / (2w - ^rx)0i0 2 , 

0lV 3 2 + / 2 V 3 0i = / 0102 - / 0102 + / (2w - ^7^)0102 • 

Proposition 25. Suppose we have an r tensorfield W0 and an r — 1 tensorfield ^0. 

( V 1 *"^ Va^Va,.^ + / V A ^W.. V V 1 -^- 1 = - / (V + r/) (1) (2) 0. 



By using the Bianchi equations and the integration by parts formula in Proposition 24 



and 25, we can prove the following energy estimates. The main new ingredient is the obser- 
vation that a can be renormalized away in the energy estimates. We can therefore control 
(3, p, a, (3, a without any knowledge of a 
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Proposition 26. The following L 2 estimates for the curvature components hold for i < 2: 

7«iTr|2 



£ / iv^i t + J2 / i v ^i 

*e{/3,p,a,/3} *6{p,ct,/3,q} ^S" 

]T /" |VN>| 7 + £ / |V**|* + esC(0o). 

*e{/3,p,tr,^} *e{p,cr,/3,a} 

+ / v^v^v* 4 * 

+ / V^V^ITV* 4 * 

J DlL.1L „" I „- I „' I »' I 1 



Proof. We first consider the following sets of Bianchi equations: 

V 3 /3 + trxp = Vp + 2w/3 +* Va + 2 X ■ P + 3(rjP +* V°), 

V 4(T + hr X a = -div *(3 + l ~x * a - ( ■* P - 2 R ■* 0, 

3 1 
V 4 p + ^tr X p = div p - -x ■ a + ( ■ p + 2-q ■ 0, 

In the V4 equations for a and p, the term a appears. This is undesirable because we do 
not have good estimates for a. Instead we renormalize both a and p as in the proof of 
Proposition 20 Define 

1 , » 
& ■■= & + -x^X, 

1 . „ 

Now, we write down the equations for a and p instead of a and p. Since V 4 x = —a + ipip, 
the equations now do not contain a. Schematically, we have 

V 4 <t = -div */3 + ipV4> + + ?p 3 

V 4 p = div p + ipVip + ip^ + ip 

where \& G {P, p, a, (3, a} is now one of the acceptable curvature components. Rewrite also 
the first equation as 

V 3 /3 = Vp +* Va + # + ipVip. (55) 



(54) 



Applying Proposition 25 yields the following: 



/ <P, V 3 /3 > 7 = / < P, Vp +* Vcr > 7 + < /3, ^ > 7 



- < div /3, p > 7 + < div <r > 7 + < /3, + ?/>W> > 7 

- < V 4 p, p > 7 - < V 4 <7, (X > 7 + < ^ + ifrVlf) + V' 3 > 7 • 
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Applying Proposition 24 then yields the energy estimates for the 0-th derivative of the 
curvature: 



r m 2 + P 2 +° 2 



< I W y + [ p 2 + t 2 + [ ^(^ + va^ + v> 3 )- 

Jh u , Jh^, Jd u ,u 

We now commute equations (54) and (55) with V* to get 
V 3 V*/3 - VV'p - VV'tf 

«l+«2+«3+i4+i5=*-l *l+i2+«3=i 



and 



and 



where 



V 4 V*p - div V*/3 

ii+i2+i3+U=i 

+ V h ^ i2 V i3 K i4 V i5 (5 + Y V^V* 3 ^), 

V 4 V i a - div *V*0 

n+*2+*3+M+«5=''-l h+i2+h=i 

(div V*/^,..^ = (V i+1 /3) B A,..A lB - 



We integrate by parts using Proposition [24] and [25] as in the 0-th derivative case to get 



DlL. U »' I I „■ I n * I 1 



21+22+^3+^4=2+1 
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We then consider the following set of Bianchi equations. Notice that a does not appear in 
any equations in this set. 

V 3 a + hrx* = -div *P + l -x •* a - ( •* §_ - 2 V •* 0, 

3 1 
V 3 P + 2 tr XP = -div 0--x-a + C-0-2ri-0, 

V 4 /3 + tr X §_ = -Vp +* V<7 + 2ouf3_ + 2x • - 3(vP ~* w), 
From these we can derive the following estimates: 

/ |V*(a>,<t)I? + / |V*£|? 
< ( |V l (p,a)|2+ I \v i §\ 2 1 + 1 V** ^ V^V^V* 4 * 

+ / V*^ V^V^KV* 4 *. 

"^"■li ii+i 2 +j3+«4=«-l 

Finally, we look at the following set of Bianchi equations: 
V 3 /3 + 2trx/3_ = -div a - 2w/3 + 77 ■ a, 

V 4 a + \trxa = -V®/3 + Aua - 3(xP ~* X?) + (C - 4r?)®/3 
Again, we notice the absence of the a terms. Hence we have the estimates: 
/ |V^ + / |VV*|? 



JH-u' Jd u ,u il+ i 2+ i 3+ i 4=i 



h +12+13 +»4 

We have thus established that 



can be bounded by the right hand side in the statement of the Proposition. To conclude, we 
use the fact that 

= if + ^ 
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and thus 



E 

*G{p,ct,/3,q} 



r 


1 \~7iiTf 1 

V w 






L 


|V^| 






L 


|V**| 







E / i v ^i?+ E iiv^v-^iu^j 



*e{/5,6-,/3,a} " — a u+i2<2 

S_ E 

*e{/5,(T,/3,a} 



□ 



Finally, we prove the boundedness of the norm 7?.: 
Proposition 27. There erases eo = e (C\),7^o) such that for every e < e 0; 

H<C(O ,1l ). 
Proof. Assume bootstrap assumption: 

K < A 2 . (56) 



Thus we can apply Propositions 16, 18 and 23 By Proposition 26, we have 

2 

E( E / ( v ^ 2+ E / ( v,lI/ ) 2 ) 

<L7^ + e5C7(O ) + / E V ^ E V^V^V* 4 * 

■/fu.u j<2 ti+t2+i3+*4<2 

+ / E V ^ E V il ip i2 V h KV ii ^ 

JT>u,u j< 2 j 1+ j 2 _|_j 3+ j 4 <i 

+ / E V ^ E v^'V^v 44 ^ 

<CK 2 + CJ2 [ (V^) 2 + E / (V'^V^V***) 2 
+ E / (V^V^ifV* 4 *) 2 + E / (V il V 42 V i3 VV i4 V) : 

H+«2+«3+M<l ^".Ji n+«2+i3+«4<3 ^ X, «>« 

It is helpful to first note that 

Eh^ii^^)^^' 



(57) 



i<2 



since V 1 ^/ can be estimated either in L 2 (H U ) or L 2 (H^ U ). By Sobolev Embedding in Propo- 
sitions [9] and 10, we thus have 

E W^WlUIl^s) + \ M\lilil°°(s) < Cei-R, 

i<l 
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A similar argument shows that 

E H^VH^a,,*) + E HVVlUsLiL*(5) + X) ii V Vlk^( S ) < ce^a 3i2 , 



i<l 



t<3 i<2 

Thus by Proposition [23} we have 

E ll v Vll^p u , a ) + E IIvvIUsl^(5) + E IIVVIU^|^(5) < c{o ,no)e^n. 



i<3 



i<2 



i<l 



Now we estimate the first nonlinear term in (57): 

il+i2+*3+«4<2 "' X, ".ii 

We can estimate, by Proposition [16} 

E ll^v'W 3 *!!^,,) 

ii,i 2 <l,i3<l 

n<2 u '- i 2 <2 u '- i 3 <2 

<e*C(O )R. 
We then estimate the term 



E 



We have, by Propositions 16 and 18 



<e5C(0 ,^o)i?. 
Finally, we estimate 



«2<1 



i 3 <2 



E 



U+i2+*3+*4<3 



By Propositions [16] and [23} 



< 



;^sup||^||^ (5 _ ) )(E su pllv i ^llL 2 (5 u , a ))(Ell v43 ^IU 2 (^)^ 



h<3 
1 



i 2 <2 



*3<3 



<eW(O ,n )R. 



Thus 

2 

E( E / (V^) 2 + E / (V*) a )<C77^ + ei(7(Oo) + c»C7(Oo,Wo)i2. 

i=0 *e{^,p,<7,/8} ^ *G{p,cr,^,Q!} — a 
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By choosing e sufficiently small we have 

n<c(o ,n ). 



This improves (56) for A2 sufficiently large. Moreover, A2 can be chosen to depend only on 



O and 1Z . Thus the choice of e depends only on O and 7Z . □ 

4.10. Propagation of Regularity. Up to this point, we have been considering spacetimes 
with smooth characteristic initial data. In the context of an impulsive gravitational wave, our 
estimates in this Section apply to spacetimes arising from the smooth data approximating 
the given singular data. For these spacetimes, we can also prove that estimates for the higher 
order derivatives of the curvature. In view of the bounds that we have already obtained, 
this follows from standard arguments. In particular, energy estimates for the higher order 



derivatives can be derived as in Section |4.9| In this case, since the equations are linear in 
the highest order derivatives, one can use the bounds from the previous Sections together 
with Gronwall's inequality to obtain the desired estimates. We omit the details and present 
final result. 

Proposition 28. Suppose, in addition to the assumptions of Theorem^ the bound 

EEE< E iivMv**ii ,J,K- 

j<J k<K i<I ^e{p,cr,f3,a} ^e{/3,p,cr,/3} 

holds initially for some I, J, K . Then, the following bounds hold for < u < u*, < u< u^, 

EEE< E iivjv$v<*n 

,J,K- 

j<J k<K i<I ^&{p,a,/3,a} ^e{/3,p,cr,/3} 

for , < e, where e can be chosen as in Theorem^ and the constant D\ JK depends only 
on the size of the initial data norm in the assumption of Theorem^ and Di jk- 

Remark 7. Notice that for approximating solutions for an impulsive gravitational wave, the 
constants in the initial norms Dj^jk depend on n. However, the initial data satisfy uniform 
estimates for each I > 2, J > with K = 0. Therefore, the corresponding uniform estimates 
will hold in the region < u < and < u< u^. 

4.11. End of Proof of Theorem |4[ Once we have the estimates for all the higher deriva- 
tives, the proof of Theorem [4] is standard. We refer, for example, to Section 6 of [28] for 
details. 

4.12. Additional Estimates for an Impulsive Gravitational Wave. In this Subsection, 
we focus our attention on the sequence of spacetimes (M. n -,9n) with characteristic initial 
data converging to that of an impulsive gravitational wave as constructed in Section |3j The 
conclusion of Theorem [5] can be applied to (-M n ,g n ). The initial data of the sequence of 
spacetimes, however, possess an additional property that a n and its angular derivatives are 
bounded uniformly in L\. In this Subsection, we prove that in the constructed spacetimes 
{M.ni9n)i the same estimates hold for all u < with < e. Moreover, we show that a n 
concentrates around u s and if a n is initially more regular uniformly in n away from u s , it is 
also uniformly more regular away from u s for < u < < e. 

We begin by proving energy estimates for a n in L 2 (H U ) and (3 n in L 2 (H_ U ) using the pair 
of Bianchi equations for (V3« n , V4/3) n ) in d8l). 
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Proposition 29. The following estimate holds: 



H u J H 



<j |V J « n |?+ / |V\0|* + I V*(a n ,/? n ) ^ V^V^V^ 

•'Sj.' •'A.,* j 1+ j 2+ j 3+ j 4= j 

+ 1 V 4 K,/3„) ^ VV 2 V* 3 KV^, 

D u,u i 1+ i 2+ i 3+ i i= i 

where here, unlike in other places, we use ^ to denote all possible curvature components, 
including a n . 

Proof. Consider the following Bianchi equations: 

V 3 a + ^trxa = V§/3 + Aua - 3{xp +* xa) + (C + 4>q)®P, 
+ 2trx/3 = div a - 2ouf3 + rja, 

Using Proposition [25] 

< a n , V 3 a n > 7 = / < a n , V <§>/?„ > 7 + < a n , ip^ > 7 

- < div a n , (3 n > 7 + < a n ,ip^! > 7 

■ < V 4 /?n, Pn > 7 + < OC n , ^ > 7 + < p n , ^ > 7 



Integrate by parts using Proposition 24 to get that for u > u', u> yf 



/ K| 7 + / \P n \*< / \Pn\y + < (OnA).^* >7 

We use the commutation formula, and note that the special structure is preserved in the 
highest order: 

V 3 V l a - V§V l /3 



714 a 



Y v^v^v* 

«l+*2+«3+*4=« 

+ Y V il ^ 2 V <3 iT 4 V i5 ^ + Y V^V* 3 ^), 
V 4 V i /3 - div V*a 

«l+«2+«3+*4=* 

+ Y V h ^ i2 V i3 K u V l5 a+ Y V^V* 3 ^), 

U+j2+«3+«4+25=«-l h+l2+i3=i 
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Performing the integration by parts as before using Propositions 24 and 25 we have 



H u J H „, 



<! \V i a n \ 2 ,+ f |V J /3 n |2+ t V J K,/3 n ) Yl V^^V^V 44 ^ 
+ J V'(a n ,/3 n ) Yl V n # 2 V i3 iTV i4 ^ 



□ 



4.12.1. Before the Impulse. We first apply Proposition 29 for the region < u < u s . Notice 
that the initial norm 



e r\\^n\\h {So , u) du<c 

i<2 J ° 



independent of n. Therefore, 
Proposition 30. 

)<C7(0 ,^o), 



2<2 — s 



independent of n. 

Proof. By Proposition [29] and all the estimates in the previous section, we have 

L 2 (Kj) 



i<2 M ^ 



<L7(0 ,^o)(l + ||V i a n || L2(H) + / || V^Hla^ W). 

Jo 

By Gronwall's inequality, we have 

E( SU P H vi ««Hi 2 (^(o lHs )) + SU P HV^nll^^j) < c(Oo,n 



i<2 U ^ e 



as claimed. □ 

4.12.2. The Impulse Region. For an impulsive gravitational wave, the initial data represented 
by the a component of curvature and its derivative are not in L 2 . By the construction of 
Section |3j the approximating sequence a n satisfies initially 

Er +2 \\V l a n \\l HSou) du<C2't 

i<2 J ^s 

In the following Proposition, we show that this bound is propagated. 
Proposition 31. 

^(sup||V i a n ||i 2( ^ feiMs+2 -„ )) + sup ||V^||| 2(S j) <C(e> ,ft )2f, 

j< 2 u ^ € u s <u<u a +2~ n 

where C(O ,1Z ) is independent ofn. 
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Proof. The proof is the same as Proposition 30 except that the initial data for V l a can only 
be bounded by 2 f in L 2 (H (u s ,u s + 2~ n )). □ 

4.12.3. After the Impulse. For u s + 2~™ < u < u^, since the initial data satisfy 



E, 

i<2 J ^s+2~ 



\V l a 



du < C, 



independent of n, we again have a better estimate. An extra challenge arises from the fact 
that the estimates derived in Proposition 31 depends on n. This can be overcome if we allow 
a loss in derivative: 



Proposition 32. 



^(sup \\V l a n \\ 2 L2{Hu{ ^ 



+2- n ,M„)) + SU P 

u,+2- n <u<u, 



^(3 n \\ 2 L2{Ku) )<c(o ,n ) 



i<l 

independent of n. 

Proof. To bound the left hand side, we can apply the argument as in Proposition |30| to get 
J](sup llV^lli^^-^)) + sup \\v l p n \\ 2 LHHJ ) < C(O ,n )2^ 

li<U* U a +2~ n <U<U tk 



i<l 



The problem in proving Proposition [32] with a bound independent of n is that according 
to Proposition 31, 1 1 V*/3 n | \l 2 (h ) m ay grow in n. However, we can still obtain the desired 
bounds for a n and Va„ (note that the argument below will not apply to V 2 a n unless extra 
assumptions on the initial data are used) by using Proposition 17, which implies that 

J2^p\\^^n\\ms^)<c(Oo,n ). 

1<1 



Taking u = u s + 2 n and integrating, we thus have 



L 2 (H 



u,+2 



-n) 



<c(o ,n ). 



i<i 



□ 



We can also prove estimates for higher derivatives away from u s . The influence of a 
propagating curvature impulse results in a loss of derivatives in such estimates compared 
to the assumed regularity of the initial data. By the remark after Proposition 28, we have 
estimates, uniform in n, for the V3 and V derivatives for ip and We now show that we can 
also have uniform in n estimates for the V4 derivatives for ip and \l/ for u>u s + 2~™. By the 
Bianchi equations Q, this in turn can be reduced to showing estimates for V4 derivatives 
of a. This is proved in the following Proposition: 

Proposition 33. For u s + 2~ n < u < e, if the initial data set satisfies 



EE E 11 

k<K i<I *e{p,o-,/3,«} 

k<K i<i - *e{/3,p,o-,^} 



v:v^„|| L2( ^ ) 



< c 



I,K- 



|V 4 fc V^„|| L2(S0i j<^. 
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then 

E E (^P\\^ i a n \\^(H u (u s +2-n,e))+ SU P 1 1 VjV\0 n | \l*(Hj) < C' I>K , 

where C^k is independent of n. 

Proof. The main step is to prove that on H u +2 -«, we have the estimate 

||V 4 fc V^|| i2( ^ +2 _„)<C7 

independent of n. For k = 0, we can commute with angular derivatives and integrate along 
the u direction using the Bianchi equation 

V 3 /3 + trx/3 = Vp + 2u(3 +* Va + 2x ■ f3_ + 3(np +* na). 

Since the right hand side has one more angular derivative, this integration loses a derivative, 
allowing us only to prove 

E HV^nlU^^n) < C 
i<I-l 

Now we consider also the V4 derivatives. Differentiate the equation by V 4 Vj and commute 
[V3, V4 Vj] on the left hand side. Moreover, except for the component a, whenever we see 

7*1 



V 4 1 \l/, we substitute an appropriate Bianchi equation from (|8j). Then we get an equation 

V 3 V^V i /3 = Vt 1 V i+2 f3 + ... 

where ... denote terms that are lower order in terms of derivatives. We then integrate this 
along the u direction. Thus 

nv 4 fc v^ n |i i2( ^ s+2 _ n) < ciivrv +2 /3 n ii i2( H as+2 _ n) + ... 

Inducting in k, we get that 

||V^V J /3 n || L 2 ( ^ +2 _ n) < C7||V l+2 ^ n || L 2 ( ^ +2 _ n) + ... 

Now have reduced to the can with only angular derivatives. By the above, we need i+2k+l < 
I. Thus we can prove 

E E \\^P n \\ LHK ^ +2 „ n) <C. (58) 

i<I-l k<mm{K,[^f±\} 

The conclusion thus follows from a standard energy estimate type argument as in Proposition 
and using (58) as the initial data on H_ u +2 -n □ 
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Combining the estimates in Propositions 30, 31, 33 and Cauchy-Schwarz, we have the 
following uniform L\ estimate. 

Proposition 34. 



SUP E / \\^ la n\\L2(S u ,u) d U < C/> 

where Ci is independent of n. 



Proposition 34 is crucial in showing that the limiting spacetime will have a defined as a 
finite measure with a singular atom at u s . 
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5. Convergence 

In this Section, we show that a sequence of initial data satisfying uniform estimates of 
Theorem [3] with converging initial data gives rise to a sequence of converging spacetimes. 
The convergence will be understood as follows: the spacetime will be identified in the system 
of double null coordinates (u,u,9 l ,9 2 ) and convergence will be established for the sequence 
of the corresponding spacetime metrics. In view of the quasilinear nature of the Einstein 
equations, we can only hope to prove convergence of our approximating spacetimes in a norm 
with one derivative fewer than the a priori estimates that we established for them. 

To get estimates of metric, we use the fact that the metric components satisfy inhomoge- 
neous transport equations with right hand side expressed as Ricci coefficients. The estimates 
for the difference of Ricci coefficients and curvature components are derived by considering 
the system of difference equations obtained from the original system of transport, elliptic and 
Bianchi equations. As was the case for the a priori estimates, the challenge for the difference 
system is a lack of any a priori control of the difference of the singular curvature components 
a. For the approximating sequence of initial data, a n is not a Cauchy sequence in L 2 (H ). 
In the proof of the a priori estimates, we handled the lack of information of the a component 
of curvature via a renormalization procedure and the observation that the system satisfied 
by the Ricci coefficients and the renormalized curvature components can be closed without 
any reference to the a component of curvature. Potentially, this property may fail when we 
consider the difference system. It is however a remarkable fact as we will show below that 
the difference equations still possess the same structure. 

The following is the main Theorem in which we estimate the difference of the metrics, 
Ricci coefficients and curvature components of two spacetimes: 

Theorem 5. Suppose we have two sets of initial data (1) and (2) satisfying the conditions 
in Theorem with the same constants C and c. By Theorem 0, we can solve for vacuum 
spacetimes (M.^\g^) and (Ai^ 2 \ g^) corresponding to the initial data sets (1) and (2) in 
the region < u < and < u < for u*,u^ < e . Let (u,u, O 1 ^ 2 ) be the coordinate 



system introduced in Section 2.2 such that the metrics take the form 

g ii) = -2{Q {i) ) 2 {du ®du + du® du) + ( 7 «) AB {d6 > A - {b^) A du) ® (dO B - (b^fdu), 

where Q = 1 and b A = on Ho and H_ Q . We can now identify the two spacetimes by 
identifying points with the same value of coordinate functions. Define g' = g^ — g( 2 \ ip' = 
.0(1) —0( 2 ) and \E' / = iffW — to be the difference of the metric, the difference of the Ricci 
coefficients and the difference of the curvature components respectively. If the data satisfy 

sup |(^)Vab(w = 0)1 + sup |(— Yi^u = 0)1 < a, 
£sup||Vy|| L 2 (Sui0) +^sup||VV / || L2(S0i „)) <a, 

ip^=X,uj i<l U i<l — 

^2 ll Vl (X' w )1k°£ 2 (So,„) - a f° r some fi xed 2 < Po < oo, 

i<l 

\\^ 2 {x^,V,V)'\\l2(h ) +sup\\V 2 (x,u,V,V)'\\lz(h ) < a, 

u 

snp\\V 2 (trx,trx)'\\ L ^s u , ) + snp\\V 2 (trx,trx)'\\L^(s ,u) < a > 
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E E li v ^'lk^o, a )+ E H v ^'lk^,o) ]<o, 

i<l y*e{/3,p,o-,/3} ^e{p,a,0,a} J 

where the angular covariant derivative and all the norms are defined with respect to the 
spacetime (1). Then the following estimates hold in {0 < u < w*} D {0 < u < u^}: 

sup|(^) 1 7abI <C'a, 

u,u Ou 

E E su pii v >ik(^)^' a ' 

u,u 



EsupHV^x,^'!!^^^) < C'a, 



i<l 



sup 1 1 V 2 (x, u, V, v)'\\l2(h u ) + sup 1 1 V 2 (x, w, ^, 5)'IU 2 (S„) + SU P 1 1 ^ 2 ( tr X> tr x)'\ \l*(S u ,u) < C ' a i 

u u u,u 

E E su Pll V ^'IU|£ 2 (^)+ E supHV^'IU^^^ J < a'a, 
»<i \*e{/3,p, CT ,^} " *e{ ft (7,fta} - / 

/or some constant C depending only on C and c and independent of a. 

By iterating Theorem [5J we can reduce its proof to the region < u < 5, < u < 5 for 
some sufficiently small 5 independent of a. Most of this Section will be devoted to a proof 
of Theorem [5] for < u < 5 and < u < 5. The proof of Theorem [5] will be carried out in 



Sections 5.3, 5.4 and 5.5 



Theorem [5] implies the following convergence result: 

Theorem 6. Suppose we have a sequence of initial data that coincide on and satisfy 
the assumptions of Theorem^ with uniform constants C and c on Ho and if . By Theorem 
[7J for every initial data in the sequence, a unique smooth solution to the vacuum Einstein 
equations (A4 n ,g n ) exists in < u < u*, < u< M t for u*,^ < e and the metric takes the 
following form in the coordinate system (w, w, ,6 2 ) : 

g n = -2(tt n ) 2 (du ®du + du® du) + ('y n ) A B(d6 A - (b n ) A du) <g> (d6 B - (b n ) B du), 

where f2 = 1 and b A = on Hq and H_ . Denote the Ricci coefficients and curvature 
components by ip n and \I/ n respectively. Identify the spacetimes in the sequence the value 
of the coordinate functions (u, u, O 1 ^ 2 ). Define also g' n = g n — g n ~\, ip' n = ip n — ip n -i and 
K :• V 

d 

sup|(^) i (7Ai?);(M = 0)| < a n , 

E E^pII^™'^ 2 ^) ^ a «> 

E ll V XX> W )nllz£°£2(S , H ) ^ a « f° r SOme fi Xed 2 < PO < OO, 
i<l 

||V 2 (x,w,r/,r/)' n || L2(Ho) + sup||V 2 (irx, irx)UlL2(s„,„) < ««, 
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li Vi <IU|i 2 (So, s )+ Yl HV i <IU^(S„,o) 
v *e{/3,p,<j,/3} *e{p,o-,/3,a} y 

/or some a n such that a n < oo, then the spacetime metrics converge uniformly to a 
continuous limiting spacetime metric g^ 

#00 = -2(fi 00 ) 2 (<iM ®du + du® du) + {^ 00 )AB{d6 A - (b^du) <g> (d6 B - (b^du) 

in the region < u < u*, < u<u^. Moreover, 

d d d d 

(ggdn, ^#«) converge to (t^#oo, g^9oo) m L™L™L 4 (S), 

(m 9 »> ^e 9n > h 9n) conver9e t0 { w 9 ^ Jam 9 - h 9oo) m L ~ L ~ L2 ^ 

(^9n, ^nY B ^(ln)AB)) converge to (^ 9oo , ^((7^ ^(Too) ab)) in m°L°°(S), 

{ Wdu 9 - ^du 9n) C0UVer9e t0 { Wdu 9 ^ ^du 9 - ] m L ~ L ^ L ^- 
As a consequence, in the limiting spacetime, 

d 2 d 2 d 2 

W 2900 ' du~m 9co ' du~ 29o ° e C ^ L ^ S ^ 

9 2 ^ oo oo 4 



dOdu dudu 

Remark 8. Notice that in the case of initial data of the impulsive gravitational wave, we 
have V*Xn converging in L^L 2 {S) for any 2 < p < oo, but not for p = oo. We thus take a 
sequence of initial data that converges in a topology which is consistent with V l x„ converging 
mLlL 2 {S). 

The convergence of the Ricci coefficients and the acceptable curvature components follows 
from Theorem [5j In Section 5/7 , we will show that the convergence of the Ricci coefficients 



and the acceptable curvature components imply the asserted convergence of the metric and 
the regularity property of the limiting spacetime. The above convergence theorem is strong 
enough to show that the limiting spacetime is a solution to the Einstein equations: 

Theorem 7. Suppose all the assumptions of Theorem [6| hold. Then the limiting spacetime 
metric satisfies the Einstein equations in L^°L^°L 2 (S'). 



This will also be proved in Section 5.7 Moreover, this limiting spacetime solution is the 



unique solution to the vacuum Einstein equations. 

Theorem 8. The solution to the Einstein equations given by Theorems [6| and [?| is unique 
among spacetimes that arise as C° limit of smooth solutions to the vacuum Einstein equa- 
tions. 
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A more precise version of the uniqueness theorem is formulated as Proposition 73 and will 
be proved in Section 5.8 Moreover, if the initial data is assumed to be more regular, the 
limit spacetime metric is more regular: 

Theorem 9. Suppose, in addition to the assumptions of Theorem^ the bounds 



EE( E 

j<J i<I #e{p,o-,/3,a} 



\viv i y n \\ L 2 {Ko) + J2 

*e{/3, P , CT ,/3} 



|V*V^ n |UW< C 



hold uniformly independent of n. Then 



E E 

j<J+2i<min{I,j-2} 



E E 

j<J+2 j<min{7j-2} 



This will be proved in Section 5.9 



For general initial data satisfying the assumptions of Theorem |3j using the construction 
in Section [3j there exists an approximating sequence of smooth initial data satisfying the 
assumption of Theorem [6] Thus, the combination of Theorems |4| [6j [7| [9] and [8] together 
imply Theorem [3} 

The remainder of this Section will be organized as follows: Theorem |5]is proved in Sections 



5.1 5.5 After the definition of the norms in Section 5.1 the proof is carried out in three steps: 



STEP 1 (Section 5.2): The difference of the metric components are estimated assuming the 
bounds for the difference of the Ricci coefficients. 



STEP 2 (Section 5.4): The difference of the Ricci coefficients are controlled assuming the 
estimates of the difference of the curvature components. This relies on the transport equa- 



tions for the difference quantities derived in Section 5.3 



STEP 3 (Section 5.5): Finally, the bounds for the difference of curvature components are 



obtained, closing all the estimates for Theorem [5j 



In Section |5.6[ additional estimates are derived for the metric components. These ad- 
ditional estimates will be used together with Theorem [5] to construct a limiting spacetime 
and to obtain Theorem [6] in Section 5.7 In this Section, Theorem [7] is also proved, showing 
that the limiting spacetime satisfies the Einstein equations. In Section |5.8[ we formulate 



and prove a precise version of Theorem [SJ establishing uniqueness of the limiting spacetime. 
In Section 5.9 Theorem [9] is proved, showing additional regularity in the spacetime with 
more regular initial data. Finally, in Section 5.10 we return to the case of an impulsive 
gravitational wave and obtain extra regularity properties for these spacetimes. 



5.1. Norms. We begin the proof of Theorem |5j Define the following L 2 norms for the 
difference of the null curvature components or their renormalized versions: 

U' = J^sup \\V l (l3,p,a,/3)'\\ L 2 {Hu) 
i<i u 



^sup||V l (p,<M, a)'\\ L 2 {Ii j 
i<i * 
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Define the following norm for difference of the the Ricci coefficients: 

°' = Yl sup 1 1 V *(*' k°£ 2 (s« A ) + Yl sup 1 1 ( V *( tr ^' V, V, w, X, ^x)'!!^^). 

These norms will be used to estimate the difference of the null Ricci coefficients except for 
those involving the highest derivatives, for which the estimates are weaker. We therefore 
also introduce the norms: 

& = sup ||V 2 (x,w,^,^)'||L2(H tl ) + sup||V 2 (x,w,r/,r/)'|| L 2 G fy + sup 1 1 V 2 (trx, tr%)1 ^(s^), 

u u u,u 

Notice that the norms 1Z', O' and & are difference counterparts of the norms TZ, O and O 
for the difference quantities. Note, however, that the former provides control of one fewer 
derivatives than the latter. This is due to the fact that convergence will be proved in a norm 
which is one derivative weaker than the corresponding norms for the a priori estimates. 
Define also the following norms: 

O" = sup(| | {x,uJ,V,V,tT X ,tTx)'\ \l*l°°{s u ,u) + \\V(x,u,V,V,trx,tex)'\\LlL*(s v , % ))- 

u 

It follows by Sobolev Embedding in Proposition [9] and 10 that 
Proposition 35. 

o" < c(6' + o'). 

5.2. Estimates for the Difference of the Metrics. In this section, we show that the 
difference of the metrics and their coordinate angular derivatives can be controlled by the 
O' and O' norms. Recall that the metrics take the form 

g® = -2{Q^) 2 {du ®du + du® du) + {^ ) ) A B{d9 A - (b^) A du) ® (d9 B - {b^) B du), 

for i = 1,2. 



Proposition 36. 

Proof. Recall that 
Hence 



supKn',^- 1 )')^)! < c8?o". 



1 1 T s~i 

u = --V 4 logfi = -nv^ 1 = -—Q^ 1 . 

2 2 2 au 

1 r5 

By integrating along the u direction, noticing that Q' = on H_ , and using Cauchy-Schwarz 
we get 

Kir 1 )'] < c5*o". 

In order to get the estimate for Q', we note that 

(o' + n^y 1 = (^ (1) r l = (IT 1 )' 2 ) + (Jr 1 )'. 

Therefore, 



□ 
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In view of the upper and lower bounds of Q in Proposition |2j 

\n'\ < cs^o". 

Using the estimates for Q', we also have estimates for 7'. 
Proposition 37. 7' satisfies the following pointwise bounds: 

sup KtU ((7-y B y)MI < a + C&0" n . 

Proof. The components of 7 solve the following ODE: 

d 

— 7ab = 2VLxab- 
ou 

This implies 

d 

— log(det 7) = fltrx- 
From this we can derive an equation for (det7)': 



^(det 7 )' = 7 - l77 ^y M ^-(det7)'^(det7) w + (det 7) w (det 7 )^(fitr X )' 

By Proposition [3j we have uniform upper and lower bounds on det 7 and uniform estimates 
for f2 and J^(det7). The previous Proposition gives I ft' I < C5-2O". Moreover, / \tix'W < 
SO" by definition. Thus, 

I (det 7)' I < C5*0". (59) 
We can also derive an equation for (jab) 1 '- 

d 

-q^Hab) = 2(Qxab) ■ 

Thus 

\(1ab)'(u,u)\ < \( lAB y{u,0)\+C r\{xAB)W + C6 1 *0"\ X AB\. 



By assumption, \{jab)'{u, 0)| < a. By Proposition [3J matrices 7^ and 7^ are uniformly 
non-degenerate matrices. Hence, regardless of whether we define the L°° norm with respect 
to 7M or 7^, we have 

\(Xab)'\ < CsupUx'IU-^)- 

A,B=1,2 u ^ 

Therefore, 

\{iab)'\ <a + C6*<D". (60) 
Now, (59) and (60) together also imply the pointwise bound in coordinates for 

|((7 _ y B )'l <a + C5-*0". 

□ 

This estimate allows us to conclude that the LP norms defined with respect to either metric 
(1) or (2) differ only by a + C5^0". 
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Proposition 38. Given any tensor, we can define its L P (S) norms \\4>\\lp(s u u) an< ^ 1 1^1 ^lp(s u 
with respect to the first and second metric. Suppose 

LP(S u ,u) - 

Then 

|W _iuii( 2 ) i ^ r>(„ j. xhrt>\\\A\\W 



In particular, is also in LP with respect to the second metric. 

Proof. This follows from the pointwise control of (■Jab)', ((7~ 1 ) AS )' an d (det7)' in coordi- 
nates. □ 

In a similar manner, we can control the (gp-) derivatives of 7'. 

Proposition 39. 

sup \\(J^W ab \\l H s^) <Ca + CSiO", 



where L 4 (S) is understood as the L A norm for a scalar function and by Proposition 38 can 
be defined with respect either to the metric y^ or y( 2 K 

Proof. By the equation 

d 

— "iab = Z^Xab 
au 

and the assumption, we have 

\\(^Yiab(^)\\l H s uaI) ^ cro + cr / a (ll(^a)*" / IU*CflU al ,) + IK^)VIU*<^)) dl*'- 

(61) 

Since 



V(logfi) = ^(v + v), 



we have 



K^cW < c\(v,v)'\ + c\r'\\n\ + c\r\\n'\ < c\( v ,vy\ + c\^i AB \ + c\n'\, 

where T is the connection coefficients on the spheres with respect to 7. Moreover, by Cauchy- 
Schwarz, 



Thus by (61), we have 

d 

\\{^QcnABK^mW{s u> u) 

<Ca + C5$0" + c£\ \if, rf, X', ^dIab, «1 W { s u>y ,)dy!. 

Using already established estimates, Gronwall's inequality and Cauchy-Schwarz, we thus 
have 

\\(^ch'AB(u,u)\\LHs^) < Ca + CS^O". 

□ 
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A consequence of the above Proposition is the estimates on the difference of the connection 
coefficients: 

Proposition 40. 

sup||r / || X 4 ( ^ ji) <C7a + C7^0". 

u,u 

Proof. This follows from the fact that T' can be expressed as a linear combination 7' and its 
first angular derivative. □ 

This implies that and are comparable in the following sense: 

Proposition 41. Let p < 4. Suppose cf) is a tensor. Then 

1 1 (V«)0 - (V< 2 ))0| \ LP{s ^j < C(a + 1 10| U^). 

Proof. Since we have estimates for •jab and (7~ 1 ) AB . It suffice to estimate 

(V (1) ) a 0-(V (2) ^ 



a 

ae A ae^ 



Using Proposition [40} we have 

ae A ae A 

^lir^-r^iiipfs^) 

<C\\V'\\ L 4( Su ,u)\\ ( t ) \\L'^(S U: u) 

<L7(a + ^O")||0||^(5_). 



□ 



We now show that the difference between the angular covariant derivatives of two tensors 
is comparable to the angular covariant derivative of the difference: 

Proposition 42. Let p < 4. Suppose <j>^ and (f)^ are tensors defined on spacetimes (1) 
and (2) respectively. Then 

||V«(0') - (V0)'|U, ( 5„, a ) < C(a + S^O")M\ l ^ Sua) . 
Proof. It is easy to see that 

v«(0') - (vcf>y = -v (1 V (2) + v (2 V {2) - 



The Proposition thus follows from Proposition 41 □ 



We now estimate the coordinate angular derivative of Q' in L A (S): 
Proposition 43. We have 

S np\\J^Q'\\ LH s^j<Ca + C5^0". 

Proof. Recall that 

1 1 X c) 

u> = --v 4 io g fi = -nv^- 1 = -j^r 1 . 

2 2 2ou 

Hence 

1 d 
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Thus, in order to estimate gp-^', we need to estimate gp-w' in L^L 4 ^). Since by the a 
priori estimates in the previous section we have ||Va;| \l°°(S) < C, it suffices by Proposition 
[42] to estimate 

/ iiW|| X 4 (s 

i<l J ° 

By the definition of the norm and Cauchy-Schwarz 

r\\&<*t\Wv^)dd< ester. 
i<i Jo 

The conclusion thus follows. □ 



The above Proposition allows us to estimate the second coordinate angular derivatives of 
i in L 2 (S). 



Proposition 44. 



rj2 

\\(^r^nhAB\\L H S u$!L) <Ca + C5?{0" + &), 



where L 2 (S) is understood as the L 2 norms for scalar functions. 
Proof. By the equation 

d 

and the a priori estimates in the previous section, we have 
d 2 

)7abKm)IIl 2 (s^) 



de c de D 



o 



The last term has already been estimated in the proof of Proposition 39 



d , d , 

) Q '\\lHs,) + \\{ 



Since 



<Ca + C6*0" 
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we have 

d 2 

d d 

^ C \\qq(v,V)\\l2(s Ui!iI ) + c \\qq T '\\l^s u ^)\M\l^(s u ^) + CI I ^ r l 1^(5^)11^1 1 U-^,) 

d d 
+ C\ \T'\ \ms u ^) 1 1 j^l \lhs u ^,) + C\ \T\ \l°°(s u ^,) 1 1 IU 2 (s Uia ,) 



i<l i<l i<l 



Using Propositions 43 and the definition of the norms O" and O', we have 

pu f) f) 

£ I (ii^M)'!!^,) + li(^)T'iU 2( 5_,) + 11(^11^,))%' 



<c5i(o"+a;)+c /"ii(^ f )7^Km)iu 2 (^,)^ / )- 



Moreover, 



Thus by ( 62 ) , we have 



;)i'ab{u,u)\\ l i {Su&) 



'd6 c d6 D ' 

32 



<L7a + C5* (O" + &) + C jT 1 1 {-Q^h'ABiu, u) 1 1^(5^)^')- 



By Gronwall's inequality, we thus have 



:^^)7^IU 2 (5„, a ) < Ca + C5\{0" + &). 



The proof of the previous Proposition also shows that 
Proposition 45. 

sup 11(^75)^11^(5^) < Ca + C^((9" + £>') 



□ 



Proposition 44 implies the following estimates of K', the difference of the Gauss curvatures: 
Proposition 46. K' satisfies the following bounds: 

sup \\K'\\ L , {Su ^ <Ca + C5kO" + &). 



Proof. This follows from the estimates of the derivatives of 7' in Propositions 39 and 44 □ 



Moreover, Proposition [44] allows us to conclude that the covariant second derivative with 
respect to connection (1) and connection (2) are comparable: 
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Proposition 47. Let p < 2. Suppose (f) is a tensor. Then 

||(v«)V - (v^^iu^) < c(a + tio") llvviU-c^)- 



i<l 



Proof. By Propositions 40 , 44 and 46 □ 



This implies that taking the covariant second derivatives of the difference is comparable 
to taking the difference of the covariant second derivatives: 

Proposition 48. Let p < 2. Suppose <f>W and are tensors defined on spacetimes (1) 
and (2) respectively. Then 

||(VW) 2 (0') - (V'miLns^) < C(a + 8*0") ^ ||V>|| i4( ^ } . 



Proof. The proof follows as that of Proposition 42, using Proposition 47 instead of Proposi- 
tion su □ 



We now prove the estimates for b'\ 
Proposition 49. 

sup 1 1 (b A )'\ \ l ~ {Sua) , sup 1 1 (b A )'\ \ Li{Sutu) , snp 1 1 ——{\A)'\\^ {M <C{&' + &). 

u,u u,u Ou u,u Ou Ou 

Proof. Recall that h satisfies the transport equation 

db A 



du 



The Proposition thus follows from differentiating the equation by J|, integrating in u and 
using the the bounds in Propositions 36, 43 and 45 and the definitions of the norms. □ 

5.3. Transport Equations for the Difference Quantities. In this Subsection, we show 
how to derive estimates for the difference 4>' of the quantities (j)^> and (j)^ satisfying transport 
equations in the V4 , V 4 or V3 , V3 directions. This should be compared with the 
previous Subsection where we considered difference quantities satisfying transport equations, 
but unlike in this Subsection, the corresponding transport equations held with respect to the 
derivatives rather than V4, V3, which themselves depend on the spacetimes under 
consideration. 

In order to estimate the quantities 0' associated to the difference between that of different 
spacetimes, we need to derive a transport equation for <f> from the transport equation for 0. 

Proposition 50. Consider spacetimes satisfy the hypothesis of Theorem^ and [3} Let <fi be 
a (0,r) S-tensor. Suppose V \4> = F. Then 

vlV + ^v 4 + n~\T l x)'<p. 

Similarly, suppose V^cp = G. Then 

\l 1 8(p£ — 
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Proof. We write in coordinates 

where Ai denotes that the original Ai in the z-th slot of the tensor is removed. This equation 
holds in both spacetimes (1) and (2). Then, we derive the equations for <ft'. We will write 
schematically without the exact constant depending only on r: 

(v 4 ) ( V 

- (n- 1 )'^ + (7- 1 ^- 1 x)> (2) 
=f' - |^(V4) (1 V (2) + (n-YH-y-W**. 

For the V3 equations, we write in coordinates 

v 3( p Al ..A r =n- 1 — ( p Al ... Ar + n^b 3 ' 



- {i- l f D ^ l x CD <P Al ... AlC ... Ar + ^ 1 b B r c BAi <p Al ... AiC ... Ar . 

We derive the difference equations as before. However, e 3 in the respective spacetimes are 
not parallel to each other and we therefore have extra terms. 

(v 3 ) (i y 

(cr 1 )^ + rjr^w&w^- - (7^) (1) (^^) (1) x {1) + (fi-W&wr^ {<P (l) - <P (2) ) 

au ad — J 

(fi -i)(D^ + (Q-i)(D 6 (i)^ _ ( 7 -i)W(n-i)W x d) + (n-^WftWr^A 

au ad — J 

«l-H*)^- + (fi-l)(2)fe(2)_| _ (7 -l ) (2) ( j ] -l ) (2) x (2) + (fi -l)(2) 6 (2) r (2)^ (2) 

au oO — J 

)' V7(l)j.(D ZO-l\(2)r / 5 



=C - Tjpryji) v«0 (1) - (^- 1 ) (2) &'^0 (2) + (n- 1 )«(7-W (2) - (n-W&r)'^. 



□ 



This, together with the commutation estimates in Proposition 12, gives the following 
estimates for (V l 0)„. 
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Proposition 51. Suppose V40 = F. Then 

V 4 (V»'~ E (V i %,?zrV i3 F)'+ E (V^MrV^V^)' 

+ (n- 1 y( E v ii (Mrvv+ E v li (r/,#v i3 xv i4 ( 

\ii+t2+i3=i U+i2+«3+*4=« 

+ (7^X)V>. 
Similarly, suppose V30 = G. T/ien 

V 3 (V»'~ E (V i %,?7) i2 V i3 G)' + J] (V^^PV^V^)' 

+ E V !i (m) !2 V !3 G+ E V ll (^,!/) i2 V i3 xV i4 0j 

\n+«2+i3=« U+i2+*3+«4=« / 

+ &'v m + ((7~ 1 x) , + &r , )v>. 

Proof. According to the already established estimates, we have \ < Q < 2. The result 



follows directly from Propositions 12 and 50 



We now use Propositions [7] and 51 to obtain estimates from a transport equation. 
Proposition 52. Suppose V 4 = F. Then 

Eiivvi 



i<l 

< 



J2\^'\\ms u , 0) + c /E 11^)11^)^' 
i<i Jo i<% 



u'<u 
1<1 



Similarly, suppose V30 = G. Then 



i<l 

< E 1 1 w i iL» ( 5b,.) + c T E 1 1 ( v j G)i 1,2^^^' 



+ Ctfi(0» + 0' + 0')(E sup(||V 4 G|| L 2 (5iiV) + |iVViU 2 (5^,)))- 

• ^ u'<u 
1<1 

Proof. We use the pointwise estimates for 7', f2' and 6' proved in the last section. 
By Propositions 42, it suffices to estimate 

since the difference can be estimated by 

E li(vW - v(0')IIl 2( 5^) < c(o" + 6')(\m\ LHSu>!i) + \\<t>\\ L ~ { s u J. 

i<l 
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We can estimate ||(V j 0) / ||l2( Suii ) using Proposition [7] and the equation in Proposition 



Recall the formula in Proposition 51 



51 



v 4 (vV)'~ E (v ii 0M?) i2 v i3 F) / + E (v^izrv^xv^y 

U+«2+*3=« U+«2+«3+M=i 

+ (fr 1 y ( E v ii (^,!z) i2 v i3 F+ E ^ h (v,v) i2 v i3 x^ u (l>) 



+ (7 _1 X)'V>. 
We estimate term by term: 

<c r(Eii^^)ii^(^))Eii( vi2i7i ) , ii^(^)^ 

^° ii<l i 2 <l 



The second term can be estimated by 



J ° i!<li 2 <2 i 3 <l 

+ c rc£\\ vii ^^x,x)'\w { s uu ,)) 

Jo 



il<l 

x (£ E I IV- (77, ^ tr X , x) I IL 3 ^(^^)) (XI IIV^IU^,))^' 

<c<E \\(^<py\\L k ms^) + csko" + a + e/)(||v0iu 2(5tti£) + ihu^)). 

i<l 

The third term is easier to estimate since by Proposition 36, (fi -1 )' can be estimated in L°°: 
f~ E IK 0-1 )' (V il (^!/) J2 V t3 F + V^Izrv^xV^) Wms^dvf 

<c?^(o" + 0'xE sup n viF n^«,«o + E SU P ll V >ll^(^v))- 

• ^1 u'<u u'<u 
i<l i<l 

The fourth term can be controlled in the same way as the third term, since (7 _1 x)' can be 
estimated in L^L 00 : 



ElK^xWII^,)^' 

<C**(0" + 6')(E SU P H V >IU 2 (^))- 



--, u'<u 
1<1 
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Putting all these estimates together using Proposition [7j we have 



LP* (S u ,u) 



i<l 

< e 1 ivy 11^(50,,) +cj2\\(v i Fy\\ LiL z i s u „) + c5||vyiU| i2( 5 U , a) 



i<l i<l 



■ u'<u 
i<~Z 



Take 6 to be sufficiently small depending only on C but not a, we can absorb the third term 
to the left hand side to get 



i<l 



i<l i<l 



+ C<^((9" + & + 0')(E sup(||V J F|| i2(5 _ ;) + HWIU^,))). 

u'<u 
?<z 

The proof for the V3 equation is analogous. □ 



5.4. Estimates for the Ricci Coefficient Difference. The results of the previous Sub- 
section can now be used to estimate the norm O' in terms of 71' . 



In this Subsection and the subsequent Section 5.5, in order to simplify notation, we will 



omit the superindices (1) and (2) in all the background quantities, for example, V3, V, f2, (3... 

Proposition 53. There exists 5 and C depending only on the a priori estimates in Theorem 
^ (but independent of a) such that 

a < Ca + C6*K' + C5*& . 

Proof. We control O' using the null structure equations. First, we consider the Ricci coeffi- 
cients tr%, x, rj. They all satisfy V3 equations of the form 

v 3 ^ = W> + * 

such that x, u do not appear in the ip terms on the right hand side; and that /3 does not 
appear as This is important because x' an d u ' cannot be controlled by CO' in the L 2 (S) 
norms and /3' cannot be controlled by C7V in L 2 (H_). In order to estimate ip using this 
equation, we need to estimate the curvature term: 



E \\(y%'\\Li»(S) < C5k0" + 6') + C&7V. 

t<l 
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and the nonlinear terms: 

£ll(v*W))1k*<*^ 

i<l 

<C5 1 HO" + 6>) + C5 l z J2 sup(||7 , V i ^V l ^|| L2( 5 u ,j 
+ C6* Yl sup \ \V n ipV i2 ip'\\ L 2 (Sul j 



u'<u 



<c^(o // + d / + o / ). 

By Proposition [52j we thus have 

Y SU P 1 1 V l (trx, X, V)'\ \ms^) 

i<l u 

< ^||V 4 (trx, % v)'\\l*(s 0v ) + sup((9: + + C) + Otfft' (63) 

i<i 

<a + (£>" + O' + £>') + Ct^ft'. 
We now consider the equation for u, which schematically looks like 

V3CU = Uljj + X^jj + Ipljj + p, 

where ip ^ x,u are the good components as above. The term ip'u or ip'x can be estimated 
since ip' can be put in an appropriate L P (S) norm and be controlled by CO'. However, the 
most difficult terms are oo'ip and xV since ui' and x' cannot be bounded in L^ t L p (S') by CO' . 
Thus, by Proposition [52l we have 



i<l 



<CJ2 II W|U 2(50is) + C7 £ HV 41 ^^, wyik^a^) + Ctffo' + C5i(0" + 0' + 0'). 
Since this holds for every u, we can integrate in L 2 in w to get 

53l|W||i2ia ( s Uia ) 
i<l 

<CV ||W|| iSL2(5ou) + / ||V*W i2 (x^)1 kz- ( s u , jdt*' + Cft' + Ctf*(<T + & + C) 



i<l 

<a + Cft' + C5 1 ? [O" + 0' + 0' 



(64) 



since we can control V l u/ by CO' after integrating along the u direction. We can estimate 
X in a similar manner as u. x satisfies 
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where ip ^ x,u. Putting Vf] in the O norm and using already obtained estimates, we have 



£liv<*i 



i<l 



<cJ2\\^x\\lhs ^) + c \\^ 1i ^ 12 x\\lW(s^) + C5H0" + a + 0'). 



As for a/, we now integrate in L 2 in u to get 



i<l 



<Z)l|V^lkw(flb,u) + C f / HV^^x'lU^Cs^j^' + C^^' + a' + O') (65) 

i<l 

<a + C6*{0" + & + 0'). 
It remains to consider tr%, rj, u. They satisfy V4 equations of the form: 

V 4 ^ = ipip + 

where ^ can be any Ricci coefficients and \I/ can be any curvature components 7^ a, a, i.e., 
all the \I/ terms can be controlled in L 2 (H) by 72.'. Moreover, since we are now integrating 
in u, all terms ip' can be estimated by CO'. By Proposition 52, we have 

J2 \\V l (tr X , V, Il-^j < E 1 1 Vi ( tr *' ^ ^'1 +CK' + C O'dyf. 

Jo 



t<l i<l 

Therefore 



£||V*(tr X)M )'i| L2(s „, a) < a + Ctfte' + C^supO" + 0' + C). 



Kl 



Using 0" < C{0' + O') and all the above estimates, we have 

O' < Ca + CSsW + CS^(d' + 0'). 
By choosing 5 sufficiently small so that C5^ < |, we have 

O' < Ca + C6*n' + C5*&. 

Notice that in particular 5 and C depend only on the a priori estimates from Theorem [4] and 
are independent of a. □ 



We now move on to estimate O' by TZ'. Recall from Propositions 19 and 20 that O was 
controlled using a combination of transport equations for G and Hodge systems. The norm 
O' can be dealt with in a similar fashion. In order to perform this scheme, we need show 
that as was with the case of the norm O where we controlled VV> from 0, the difference ip' 
satisfies elliptic equations with 0' as a source. This is given by the following Proposition: 

Proposition 54. Let <ft be a (0,r)-tensorfield. Suppose 

div (f> = f, curl cf) = g, tr (p — h 
on each of the spacetimes. Then 

div (f)' ~ /' + (7 _1 )'V0 + V<f>, 
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curl (f)' ~ g' + <y )'V0, 
£r <j>' ~h' + (7 _1 )'0. 

Proof. This is a straightforward coordinate computation. □ 
Based on the estimates we have on 7', (7 -1 )' and T', and using the estimates in Proposition 



14 we have 



Proposition 55. Let and cf)^ be totally symmetric r + 1 covariant tensorfields on 2- 
spheres (§ 2 ,7^^) ; (S 2 ,^ 2 ^) respectively satisfying 

div = /<*>, curl 0« = g®, tr$« = h® , 

for i = 1,2. Then 

||VV|| L2( s) <Ca^(l|V7ll^(5) + l|V^|| L2( 5) + \\V l h\\ L 2 {s) + M\ L 2 {S) ) 



i<l 



+ Cj2(W^f'\\LHs) + IIVVII^(S) + HV^'ii^s) + MWlhs)) 



Proof. This is a direct consequence of Propositions 54 and 14 □ 



We can now begin to estimate O' . First, we define B' and derive estimates for 0' using 
transport equations. Let 9' denote (Vtrx)', (Vtrx)', We have the following 

estimates: 

Proposition 56. 

1 1 V^VtTX Vtrx, fJ>, ri'\\L>(s UA ) <Ca + C6* {& + W). 

i<\ 



Proof. From the proof of Proposition 20, we know that each of Vtrx, Vtrx, p, /1 satisfies 
either 

V 3 9 = TpVlpH + ^ + Ip^H, 

or 

v 4 e = fpvtpn + ^ + ^h, 



where as in Section 4.8 , we use the notation tpn e {trx, X, t], 77, u, trx}, ipn_ € {trx, t], rj, u, trx, x}, 
V H e {P,p,a,P} anHWff e {p,a,P,a}. 

We focus on those 9's satisfying the V4 equation. The other case can be treated analo- 
gously. By Proposition [521 we need to control 



1 



i<l 



i<l 

and 



L 2 (S u ,u)- 



i<l 

We first estimate the ' terms. Among those we first look at the term with two derivatives 



on the Ricci coefficient. We have, using Proposition 53 



\\{^ 2 M'\\lIl^) < ^ + L7||^|| L i^(5_) + ^i||V 2 ^|U S L 2( 5_) < C{a + & + 6fa'). 
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The other term with Ricci coefficients can be estimated by O' and thus by Proposition 53 

*1+«2<1 

We then move to the term with the curvature components: 



Y \\(^ i2 ^ i3 ^Hy\\L k L H s UA ) 

<C(a + 8*-R! + 8$&) + C8 1 * £ HW* 



i<l 



<C(a + 5»H' + 8*6'). 
The terms without ' can be estimated using the a priori estimates in Theorem |4j 



< c. 



i<l 



Thus, by Proposition [52j we have 

^HV^Vt^Vtr^/i^yil^^) < Ca + C6*(& + K'). 



i<l 



□ 



We now consider the difference quantities k' and k/. As in Proposition 22 they will only 
satisfy L 2 (H U ) and L 2 (H_ U ) estimates. As in the proof of Proposition 22, in what follows, we 
allow ip a l so to be uy and uj} . 



Proposition 57. 



i<l 



^||VV|| l2( hj < CS^a + O' + K'). 



Proof. Recall equations (49) and (50): 



or 



We focus on k. As k is easier. The only terms that are new compared to the proof of 



Proposition 56 



Ell v W)'lk^, H ) and Ell^^'lk^ 2 



i<l 



i<2 



Thus using the estimates in the proof of Proposition 56, we have 



i<l 



i<l 



i<2 



Integrating over u in L , we get 
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EIIVk'i 



L 2 (H U ) 



<C6*a + CS{& + K') + C5kJ2\^^'Hws^j + II^OM'!! W 



i<l 



i<2 



Thus, using the definition of the norms O' and 1Z', 



J3l|VV|| L2(Hu) < C5h + C5(d' + TZ'). 



i<l 



Similarly 



^||VV|| i2(ifu) < C5*a + C8{& + 11'). 



□ 



We use Propositions 56 and 57 and elliptic estimates to derive all the estimates for V 2 ip'. 
Proposition 58. 



O' <Ca + CTZ'. 



Proof. We will first prove 

sup ||V 2 (x, u,V,v)'\\l*(h u ) + sup ||V 2 (x, uj.,v,v)'\\^iH u ) <Ca + C6*d' + CTZ' 

u u 

We have the following div-curl systems: 

div (x, w, Tf, rf) = (Vtrx, p, ft) + ifnj) + {fi, p, a), 
curl (x, u, rj, rj) = (Vtrx, p, k) + 4>4> + (J), p, a), 



and 



div (x, w, T], V) = (Vtrx, fx, k) + ^ + (p, £), 



curl (x, w, 77, rj) = (Vtrx, //,«) + ^ j)j) + (p, &, 0). 

By Proposition |55j we need to control the terms involving 

SU P ^2 1 1 V *( Vtr ^' ^ K )'\ \l 2 (h u ) + sup ^ 1 1 v *( Vtl X P, «)'| U 2 (HJ) 



the curvature terms 



sup 1 1 V l (/3, p, | L 2 (Hu) + sup 1 1 V l (p, &, (3)'\\ L 2 {K j 



i<l 



and the lower order terms involving the Ricci coefficients 

X;nv*w)ii^)+E E iiv i (^v)iil 2( hj. 

»<i V>^x,w 



i<l 



The terms involving are controlled by in the Proposition 56 by 

C5*(a + & + K'). 



56 



we have 
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The curvature terms can be estimated using the definition of TZ' by 

CTZ'. 

The remaining terms can be bounded by 

x;nvw)iu^) 

i<l 

i<l i<l 

<Ca + C6*& + C6SK'. 
Thus we have 

sup \ \V 2 (x,uj,V,v)'\\l2(h u ) + sup||V 2 (x, u } 7i,v)'\\Li(H ub ) <Ca + CS^O' + CTZ'. 

u u 

Therefore, together with the estimates for V 2 tr^' and V 2 tr\;' in Proposition 

6 < Ca + C6*& + CTZ'. 
By choosing 5 sufficiently small depending on C so that C5^ } we have 

6 < Ca + CTZ'. 

Notice that C and 5 depend only on the a priori estimates in Theorem [4] and are both 
independent of a. □ 

5.5. Estimates for the Curvature Difference. In order to finish the proof of Theorem 
[5j we need to estimate TZ' by Ca. This will be proved using an energy-type estimate. 

Proposition 59. 

E( E / ( v ^') 2 + E / ( v ^ 2 

<Ca 2 + C5(n') 2 + C\\(m'\\l lLiL2{s) 
+ C E IKV' il V i2 VV i ^) / ||| SI/ii2(s) + E \\^ 12 ^ 1 ^)' wsy 

Proof. The energy estimates for the curvature difference follows from the difference Bianchi 
equations. Once these are derived, the energy estimates will be obtained as in Section [4] 
by integration by parts in appropriate subsets of the difference Bianchi system. In the 
proof below, we concentrate on the most difficult case which involves renormalized curvature 
components. The estimates for the other components can be derived in a similar fashion. 
Recall that 

V 3 /3 + trxP = Vp + 2u(3 +* Vcr + 2 X • §_ + 3(r/p +* rja), 
V 4(T + htrxa = -div *(3 + l -x * " - C •* $ ~ H ■* 0, 
V 4 p + \tr X p = div $ - -% ■ a + C ■ + 2r/ • f3, 
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As before, we renormalize the equations. Define 

1 A A 



1 

Schematically, we have 



V 3 /3 = Vp +* V<j + vpVip + + v 3 , 
V 4 <7 = -div */3 + ipVip + ipty + ip 3 , 
V 4 p = div @ + ipVip + ipty + ip 3 . 



From these we can derive equations for /?', a' and p'. By Proposition 50 



fo-iy 

V 3 /3' = (V 3 /3)' + + fi _1 (6 A )'V a /3 + n-^T" 1 ^)^ + fi _1 6r'/3, 



IV 

— 1 /_,— 1„ ,\* z- 



V 4 <r = (V 4 cr)' + V 4 a + fi- 1 ( 7 - 1 x ) / a, 



iv 



v 4 p = (v 4 p)' + ^rf v 4 p + n- 1 ( 7 - i x) P- 

Moreover, we can write schematically 

(V(p,a))' = v(p', ( r') + r'(p,a) 
(div /?)' = div + yv/3 + r'/3. 



Therefore, 



V 3 /3' = Vp +* W + + ^ + ^ 3 )' + 

V 4 a' = -div + (^V^ + ^ + + Ft, 
V 4 p' = div p + (V>VV> + + ^ 3 )' + 



where Fg, Fj and F^ satisfy the bound 

^||V i F|U2 iiL2(5 )<C^(a + ^ / ). 



j<i 



Therefore, multiplying the first equation by /3', integrating by parts and using the 
equations, we have 



<Ca 2 + 1 1 (/?', p', a') + ^ + ^y\\LiLW(s) + 1 1 03', p', <x')F AM | kzi^(S) 

<CV + 1 1 (^,^011^^(5) + ||(W + ^ + ^y\\ 2 Ll ^ LHs) + +II^IIW 2(S) 
<Ca 2 + C5(K') 2 + ||(V»VV + V* + ^yilliLimsy 
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We now look at the equation the first covariant angular derivative for the curvature compo- 
nents: 

V 3 V/3 = V 2 p +* V 2 <t + K(p, a) + + + V>W + ^V 2 ^ + </> 2 W> + 

V 4 Va = -div * V/3 + + ^ 2 W + + ^V* + ipV 2 ip + + 
V 4 Vp = div V/3 + Kp + ?/> 2 W + + + ipV 2 ip + ip 2 Vip + ip 4 . 

A similar argument as before involving writing the equations for V(/3',/5', &') and integrating 
by parts yield 

|V/3f+ I (|Vp'| 2 + |W| 2 ) 

<Ca 2 + C5(n') 2 

+ 1 1 (If* + t/; 2 V* + V^* + + VV 2 ^ + (VVO 2 + ^ 2 W> + ^)'\\l2 LlL 2 {s y 

□ 

From this one can conclude 
Proposition 60. 

W < Ca. 

Proof. By the previous Proposition, we need to estimate 

\\(m'\\i iLi L HS) , 
J2 ii(^v*w 3 *)iiw 2( 5), 

«1+12+«3<1 



w^H^miiLiv 



By Proposition 46 and the definition of the norm 1Z' and the a priori estimates from Theorem 
[IJ we have 

11(^)11^^(5) < WK'nliLims) + \\K&\\kwW < C5 2 {TZ') 2 . 



To estimate the remaining terms, we will repeatedly invoke Proposition 42 to exchange (V*)' 
and V(\I // ) (or (Vip)' and V(^/>') etc.), at the expense of an error term ?r(lZ') 2 . 

\\& h vHv i3 miiLiL* { s) 

n+«2+*3<i 

<^ 2 (^0 2 + c(^^||v^ii^ i ^ ( 5)) 2 (Eii v ^ / liW 2 w) 

ii<li 2 <2 i 3 <l 
+ C C^2 W^ 11 ^Wl™L^L°°(S))(^2 1 1 ^ I WlLlLHS) ) (/, l|V 1 '0'lli«»i«'L4(S)) 

il<l «2<1 *3<1 

<C5(a 2 + (ft') 2 ), 
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where the second term is estimated by Theorem [4] and the definition of 1Z'; and the last term 
is estimated by Theorem [4] and Proposition 53 Finally, we estimate using Theorem [4] and 
Proposition [53] 

||(V il v fe vv fe v) , ii! s ^ ( 5) 
<c5 2 (n>) 2 + c(£ liv^ll W-(^)) 2 (E IIVVII WW 



ii<l i 2 <2 



«3<2 



<L7(5(a 2 + (^') 2 ), 
Putting these together, we have 

(ll') 2 < Ca 2 + C5(TZ') 2 . 
The conclusion follows by choosing 5 sufficiently small depending only on C. 



□ 



This implies all the quantities that we have estimated in this section can be estimated by 
Ca: 

Proposition 61. 

d d 2 , ~ , 

\g\ U« ls°l-(5) , 1 1 qq9\ U» l^l*(s), I \-^9 , \\l«'L^lhs),0',0', K < Ca. 



Proof. This is a direct consequence of Propositions 36, 37 39, 43 44 45 49 53 56 58 and 
[601 □ 

This concludes the proof of Theorem [5} 

5.6. Additional Estimates on the Difference of the Metrics. In this Subsection, we 
derive some additional estimates on the difference of the metrics. They follow directly from 



the bounds already derived. This will be used in Section 5/7 to show that the limiting 
spacetime metric is in the space asserted in Theorem [6] and satisfies the Einstein equations. 



Proposition 62. 



I — 
du 



(y'M,b')\\ L * {SuA) <Ca. 



Proof. 7 satisfies a transport equation in the L direction: 

£ L 1 = 2%. 
From this we can derive an equation for 7': 



d_ 
du 



7 



-b 



1 9 y - (b A y d 



36 A 



de A 



Estimating the right hand side using Proposition |6TJ we have 

,, d „, 



du 



7 IIl 4 (s u , m ) 



< Ca. 



Q also satisfies a transport equation in the L direction: 
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As for 7', we can derive an equation for By direct estimates using Proposition 

we have 



61 



\\^'\\ LHSuA) <Ca. 

Finally, we move to §^b' . b does not satisfy a transport equation. We therefore resort to the 
equation 

f)h A 



Applying we get 



du 

du ' 



f) 2 h A d 

dudu du 



Since {b') A = on H Q , f- = on H . Thus 



By Proposition [6TJ the right hand side < Ca. □ 
Proposition 63. 

d 2 

(7 MM\\ms UA )<Ca. 



du 2 

Proof. We prove this Proposition by taking an extra derivative of the equations in the 
proof of Proposition 62 By 

^7 = 2%, 

we have 

d 2 „ d ,„ s d ,, A d 



r7 = 2-(0 X ) - ^--j). 



du 2 du — du d6 A 
Notice that 

d 

tt^-X = v 3 x - XX- 
du— — — 

Thus using the null structure equation for V3X, and estimating directly, we get 

1. a 2 .,, 



To estimate J^^, we use 



n du 2l IIl 2 (5^) < Ca. 
1/ 9 lA d 1 



Differentiating with respect to using the estimates for J^fr 4 and the null structure equa- 
tion for V3CU, we get 

1 1 d 2 , 



9m 



IU 2 (s u , 2 ) < Ca. 



Finally, differentiating 

aw 

twice in the u direction, we get 

Q3 b A Q2 



2/-A\ 



dudu 2 du 2 
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Using the null structure equations as well as the Bianchi equation for V3/3, we see that all 
the terms on the right hand side can be bounded in L 2 (S U: u). Thus we can integrate to get 



du 1 



\lhs u ,u) < Ca. 

□ 

It is also possible to prove estimates for the mixed second derivatives 
Proposition 64. 

" °~ {i,a,V)\\*{s^><Ca. 



1 dud6 A 



Proof. This can be proved in a similar way as Proposition 62, taking 7^4- instead of 
derivative. □ 

On the other hand, u derivatives can only be taken once and can only be estimated after 
taking the norm: 

Proposition 65. 

d 

||^(7,^,Mlk°L~(s Uj „) <Ca. 
Proof. We can directly estimate the right hand side of the equations 



□ 



The mixed u and 9 derivatives can also be estimates: 
Proposition 66. 

/ rV 7 /\ 



1 1 dud6 A 

Proof. Take ^p- derivatives of the equations 



d 2 

(7,fi ,b)\\ L w L4(Su ^ < Ca. 



^-(b A y = -4(f!V)', 



and apply the estimates in Proposition 61 □ 



Similarly for the mixed u and u derivatives: 
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Proposition 67. 

d 2 

\\-^(j,n,b)\\ L Po L2{Su ^<Ca. 
Proof. Take derivatives of the equations 

^-(b A y = -4(f!V)', 

and use the null structure equations. □ 

We do not have estimates for the second u derivatives of the metric, since they can only 
be understood as measures even for the initial data. However, the following substitute is 
useful in showing that the limit spacetime that we construct satisfies the Einstein equations, 
in particular, Ry^ = 0. 

Proposition 68. 

ll| I ((7- 1 ) AB | I 7^)'ll^ L . ( ^ ) <Ca. 



Proof. Notice that in coordinates, 



trx =( 7 - 1 )^|- 7AB . 



Thus we want to estimate J^tr%'. Using the null structure equation 

V 4 trx + ^(trx) 2 = ~\x\ 2 - 2wtrx, 

we have 

|-tr X ' = (-v](trx) 2 - V\x\ 2 - 2nutx X y. 
The conclusion follows from estimating the right hand side directly using Proposition [611 □ 



5.7. The Limiting Spacetime Metric. In this Subsection, we prove Theorems [6] and [7j 
We therefore assume that the conditions of Theorem [6] hold. Consider the sequence of metric 

g n = -2(tt n ) 2 (du ®du + du® du) + ( 7 „)ab(g$ A - (b n ) A du) <g> (d6 B - (b n ) B du). 



By Proposition 61, 7n , b n and Cl n converge uniformly to their limiting values 700, boc and 
Qoo and are therefore continuous functions of (u,u,6 l ,6 2 ) and define a continuous limiting 
spacetime metric 

5-00 = -2(fi 00 ) 2 (rfw ®du + du® du) + ( 7o cW(d6» A - (b^du) <g> (d6 B - (6 00 ) B du). (66) 

To complete the proof of Theorem |6j we need to demonstrate the desired regularity of the 
limiting spacetime. This follows easily from the estimates in the previous Subsections. 
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Proposition 69. 
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,d_ 
"36 



d 2 



>9n 



d 2 



Oi, 

d 2 



9ni Qy^ n ' 



d d 

converge to (—g^, —g^) in L™L™L 4 (S), 



!Jn- —9n) converge to (^<?oo, 



' dudO*'" du 2 
d 



d d 

^9n^((l n L ) AB Q^(ln)AB)) converge to 



. d 2 

^d6du 9n ' 



d 2 



dudu 



g n ) converge to 



d 2 
d 

d 2 



°~ ffoo.^oo) inL?L?L 2 (S), 



du 2 ' 



O 9<Xn rv 

du du 



9 ((^^l^ooW)) *n L~L%L°°(S), 



-9 



d 2 



dOdu ' dudu 



9oo ) m L~I%L\S). 



Proof. That 



d d 

d~e 9ni du 9n converge in L u L u lA ( s ) 



follows from Propositions 61 and 62 respectively. 
That 

£)2 f\2 £)2 

( ~ ' ' ' -g n converge in L C ^L C ^L 2 (S) 



J9n 



9r, 



de 2 " 1 " du 2 "'" dude 

follows from Propositions 61, 63 and 64 respectively. 
That 

° 9 <hn 1 ) AB 4-(ln)AB) converge in L^I^L^S) 



o 9ni rs • 

du du 



du 



follows from Propositions 65 and 68 respectively. 
That 



d 2 



-Qn, 



d 2 



dOdu ' dudu 



g n converge in L™L P U °L\S) 



follows from Propositions 64 and 5^3 respectively. 
Proposition 70. 



□ 



d^_ 

de 2 



d 2 d 2 ^ 

]9oo: 



9oo, a _. an 9oo, T^jfi'oo G C^C^L 2 (S), 



dudO du 2 ' 



9 9 ^) AB ^oo)A B )eL-L-L- 



n 9oO ) r-i 

du du 



OS), 



d 2 



-9 



d 2 



dOdu dudu 



9oo G L™L™L\S). 



Proof. To prove the first two statements, notice that since g n are smooth, the convergence 
in L™L™L P (S) implies that the limit is in C°C^L P (S). 

For the latter two statements, we use the fact that for p Q < oo, if f n — > f in L Po and f n is 
uniformly bounded in L°°, then / G L°°. □ 



Proposition [70] allows us to conclude that all first derivatives of g^ can be defined and 
belong to an appropriate space. These norms also allow us to conclude that the product of 
any two first derivatives of the metric belongs to L 2 (S) . For the second derivatives, 

Proposition 70 shows that all second derivatives of the metric except ^{j, b, Cl) are defined 

as functions in L°° L 2 (S) . The second derivative J^(7, b, fl) is merely a distribution. As 
we will see below, its definition is not necessary to make sense of the Einstein equations. 
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We now show that all components of the curvature for the limiting spacetime, except for 
RuAuB, RuuuA and Ruuuu, can be defined at worst as functions in L 2 (S) . For this we 

use the coordinate definition of the Riemann curvature tensor: 

Proposition 71. In the limiting spacetime, all components of the Riemann curvature tensor, 



defined by the coordinate expression (61) expect for RuAuB, RuuuA and Ruuuu, can be defined 
as functions in L^L^L 2 ^) . Moreover, R(L, L, L,ca) and R(L, L, L, L) can be defined as 
functions in L? (S) . 



Proof. Using (67), it is easy to see that if 5,a,n,u ^ u, the expression has at most one 
u derivative. Thus, by the regularity properties of g^, these components of the Riemann 
curvature tensor are well-defined as L^L^L 2 ^) functions. 

In the case that exactly one of 5, a, fi, v is equal to u, we can assume, by the symmetry 
properties of the Riemann curvature tensor, that 8 — u. Then, the formula for the component 
of the Riemann curvature tensor has at most one u derivative and by the regularity properties, 
it is well-defined as L 2 (S) functions. Notice, however, that strictly speaking, at the 

level of regularity that we have, the symmetry properties of the Riemann curvature tensor 
may not hold. Consider, for example, the component R\2ui- Using the formula above, we 
would have a term giA(g A -gi2,u),u, which is not defined. However, the sum of all terms 
containing two u derivatives is equal to the expression 

d 

which can be rearranged, up to terms with at most one u derivative, in the form 

d 

-^\gi a g- a gi2,u)- 

This expression vanishes since gi a g- a = 0. The above calculation is a reflection of the sym- 
metry properties of the Riemann curvature tensor and provides an appropriate distributional 
definition of RaSuu for which the symmetry properties formally hold. 

We have therefore shown that all components of curvature with the exception of RuAuB, 
RuuuA and Ruuuu are defined as L^°L^°L 2 (S') functions. However, replacing the double null 
coordinate system {u^u^O 1 ,6 2 ) by the null frame (L, L, e 1 , e 2 ), we can show that additional 
components of curvature are defined as functions in L^L^L 2 ^). To see this, we can write 

R(L, L, L, 6a) = RuuuA + b B RuBuAi 

and 

R(L, L, L, L) = R 

UUUU 

+ b RuuuB + b RuAuu + b b RuAuB — Ruuuu + 26 RuAuu + b b RuAuB- 



Note that 



Ruuuu = ~guu,uu + good terms = -^b A b B g AB ,uu + good terms, 
RuuuA = -^gAu,uu + good terms = ^b B g AB ,uu + good terms, 
RuBuA = ~^gAB,uu + good terms, 
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where good terms denote terms with at most one u derivative. Therefore, the terms with 
two u derivatives cancel in R(L,L,L,A) and R(L,L,L,L) and they can be defined as 
L™L™L 2 (S) functions. □ 

We now verify that the limiting spacetime that we have constructed satisfies the vac- 
uum Einstein equations in the sense that relative to the system of double null coordinates 

(u,u,e\e 2 ), 

iV = in L™L™L 2 (S). 

For this, we use the coordinate definition of the Ricci curvature: 

f) f) 
e> — rp rp _i_ r p r A — v p r A 

dxp Mi/ dx u PIJL px pv px pv ' 

Proposition 72. The limiting spacetime satisfies the vacuum Einstein equations in the sense 
that relative to the system of double null coordinates (u, u, 9 l , 9 2 ), 

R^ = mL™L™L 2 (S). 



Proof. Now, notice that with the metric given by (66), we have 

(g£) m = (g^r = (g^T 1 = (g^T 2 = o. 

Therefore, T^ u does not contain a term with g UUtU . In other words, the expressions for 

R uu , R u a, Rab for A, B — 1, 2 
do not contain two u derivatives. Therefore, we have 

Ruu = R u a = Rab = 

as functions in L 2 (S) . For the component R U u, there are terms involving two u 

derivatives of the metric. In particular, we have 

d d 

Ruu = tt~ TiL — 7— Tt. + terms involving at most one u derivative of the metric. 

ou - ou - 

From this expression, we see that the terms involving two u derivatives of the metric cancel. 
Thus we also have 

Ruu 

as functions in L 2 (S) . Similarly, for Rau, we have 

d_ r u _ d_ 
du A ^ du 

As for R U u, because of the cancellation, we have 

RuA = 

as functions in L™L™L 2 (S). 

It now remains to study R uu . Indeed this term involves two u derivatives of the metric 

and there are no cancellations to remove this term. We have 

d d d d 

Ruu = 7^-FL — Tq - F^ u +7^r„„+— I^ u +terms involving at most one u derivative of the metric. 
ou — ou — ou A — ou - 

The first two terms cancel each other. Note that 



Rau = ~T^ n — — T^ A + terms involving at most one u derivative of the metric. 
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since 

9uu = 9ub = 0, g Au = 0. 

We are thus left with 

It might seem that this term behaves like ^7 which we cannot control. However, this 
particular combination of derivatives, by Proposition 68, is in fact in L™L™L A (S). Thus 



R mi = in L™L™L 2 (S). □ 

5.8. Uniqueness. In the Subsection, we show that the spacetime constructed in Theorem 
[6] is the unique spacetime among the class of continuous spacetimes which arise as C° limits 
of smooth spacetime solutions to the vacuum Einstein equations. In view of Theorem [5j 
uniqueness follows easily if the spacetime solution is assumed to be in the class of spacetimes 
satisfying the conclusions of [5] Nevertheless, we prove uniqueness within a larger class of 
spacetimes for which it is not a priori assumed that the spacetime is Lipschitz, or that it is 
more regular in the u and the angular directions. 

The following is a precise formulation of the uniqueness theorem (Theorem [8]) and is the 
main result in this Subsection: 

Proposition 73. Let (Ai^\ g^) and (Ai^ 2 \ g^) be two C° Lorentzian spacetimes in double 
null coordinates, defined in < u < u*,Q < u< and taking the form 

g (i) = -2(Q^) 2 (du ®du + du® du) + (7 W )ab(^ A - (b {{) ) A du) g> (d0 B - {b^) B du) 

such that 

• Q = 1 and b A = on H and H and 7^|h = 7 (2 ^|// and 1^\h = ^liio 

• for i = 1,2, there exists a sequence of smooth spacetimes (Ain\ g^n) such that for 
every n, gn is defined in < u < and < u < and takes the following form 
in double null coordinates: 

g® = -2{ttf) 2 {du ®du + du® du) + (^) AB (d9 A - {bf) A du) ® {d6 B - {b ( $) B du) 
with the property that 
(1) 

g {l) in C° in < u < u*,0 < u< u*, 

(2) The initial data for g$ satisfy the assumptions of Theorem^ uniformly, 

(3) The initial data for g$ converges to the initial data for <?W , i.e., 

9n\H -> g {i) \h and g® \ Ko -)• \h 
in the norms in the assumptions of Theorem^ 
Then, if u^,u^ < e, where e is as given in Theorem^ 

g W = g m m 0<tt<M„0<«<M, 

Proof. By Property (3) in the assumptions, for every i, there exists rij such that 

^P \(^y(lAB)' n M = °)l + SUP ||( j^(7AflU« = 0)1 < 2~\ 

E E^Pl^VnJI^^ + ^supllVV^JI^^)) < T\ 

tpy^x.io i<l i<l — 
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E 1 1 U Yn t I k°L2(S ,„) ^ 2 ~ l fOT SOIIle fixed 2 < PO < OO, 

»<1 

u 

sup||V 2 (trx,trx)l.|U2 (5uo) + sup 1 1 V 2 (trx, tr*)^ | \l*(s ,u) < 2 ~ ? > 

u u 

E Eh v %ii^ 2 (^ h )< 2 ^ 

*G{^,p,cr,^} i<l 

By Theorem [5j there exists C depending only on the uniform bounds in the assumption (2) 
such that 

\gS-9%\<C2-\ 

Therefore, in C° norm, 

\9 {1) -9 {2) \<C\gU-g$\<C2-\ 
Since this holds for every i, we have 



g w = g {2) . 



□ 



5.9. Propagation of Regularity. Using Proposition 28, we show that if the initial data is 
more regular in the V3 and V directions, then so is the limiting spacetime. 

Proposition 74. Suppose, in addition to the assumptions of Theorem [5[ we have the fol- 
lowing estimates for the initial data: 

££( £ 11^*11 E llv^v^|| iW )<c. 

j<J i<I 9e{p,a,P,a} ^e{/3,p,cr,/3} 

Then 

EE( E su P ||viv^iu 2(iU + supiiv j 3 v^iu 2(Htt ))<c". 

j<j i<i *e{p,o-,^,a} - ^e{p,p,<T,i3} u 



Proof. By Proposition 28 , for an approximating sequence of spacetimes, the following bound 
holds independent of n: 

EE( E Bup|ivjv** B || tf 

j<J i<J *e{p,o-,/3,«} - Ve{P,p,cr,fS} " 

The conclusion thus follows. □ 
Using the equations 

1-7 = 2^X, 
= -4fi 2 ^, 

aw 



together with Proposition 74, we can show that g is in the desired space as indicated in 
Theorem M The details are straightforward and will be omitted. 
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5.10. The Limiting Spacetime for an Impulsive Gravitational Wave. In the Sub- 
section, we show that for the case of an impulsive gravitational wave, the spacetime is more 
regular than a general spacetime with merely bounded x, thus concluding the proof of The- 
orem [TJ In particular, we show that is well defined as a measure which is the weak limit 
of a n . All other curvature components can be defined in L°°. Moreover, the singularity of 
ctoo is supported on a null hypersurface and the constructed spacetime is smooth away from 
this null hypersurface. 

5.10.1. Limit of the Curvature Component a. We first show that a n has a well defined limit 
as a measure. 

Proposition 75. In the case of an impulsive gravitational wave, for every u and every 
d = (O 1 ^ 2 ) £ S 2 , aoo(u,$) is a measure on [0,u^\, defined as the limit of a n (u,"d) 



Proof. By Proposition 34, for every u and a n has uniformly bounded L\ norms. Thus to 
show convergence of a n , we only need to show that for every u and every d = (6? 1 , # 2 ) e S 2 , 

a n (u, v!, {})dyf 

converges to a limit $)([0, u)) as n — )■ oo for all u such that a^u, "#)([0, u)) is contin- 

uous. By the equation 

V 4 X + trxx = -2wx - a, 

we have 

a n (u,ri,$)duf = J (VL~ 1 —x n + tYXnXn + Zu n Xn)(u,ri,$)dy!. 
Integrating by parts and using 

d 

we derive 

a n {u,$)du = {tt~ 1 x n ){u,u,$)-{tt- 1 x n ){u,u = 0,$)+ / (tTx n Xn)(u, $)dy!. (68) 



By Proposition 61 



(tiXnXn)(u,u',$)dlJ? 

converges. Thus, in order to show weak convergence of a n , we need to show pointwise 
convergence for Vt~ l x n for all u such that ctoofu, "#)([0, u)) is continuous. By the construction 
of the data for an impulsive gravitational wave in Section [3j we know that for the initial 
data 

(Xu)ab(u = 0, u, 0) < C2~ n for u < u s or u > u s + 2~ n . 
Using the equation 

V 3 X + ^trxx = V®r] + 2ux - ^trxx + V®V, 

one sees that 

(XnW«, u, #) < C2- n for «<M s or«>« s + 2~ n . 
Therefore, Xn converges for all u ^ u s . By Proposition^} Vl^Xn also converges for all u ^ u s . 
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It remains to show that the limit a 00 (M, $)([0, u)) is discontinuous at u = u s . In view of 



(68), it suffices to show that x has a jump discontinuity across H_ u . This follows from the 
equation 

Rewriting in coordinates, we have 

d d 1 

Consider the coordinate system (u, u, e 1 ^ 2 ) such that 

u = u- b A 9 A . 

Since \b&\ < Ce, this is well-defined in a local coordinate patch. Thus, 

d 1 ~ 1 ~ 

- l" XX + 2 tr XX) = VOr? + 2wx - -trxx + V^V- 

Notice that x is continuous and the right hand side of the equation is also continuous. Since 
for the initial data, x(uo,u,9) has a jump discontinuity for u = u s , x a l so has a jump 
discontinuity across u = u s . □ 

5.10.2. Control of the Curvature Components (3, p,o~, {3, a. Except for the curvature com- 
ponent a which can only be defined as a measure, all other curvature components can be 
defined as L™L™L°°(S) functions: 

Proposition 76. 0, p,o~, (3, a can be defined as functions in L°° (S) . 

Proof. That p,a,/3,a can be defined as functions as L°° (S) follows from Proposition 



74 



and the Sobolev Embedding Theorem. To show that f3 can be defined as a L°° (S) 

function, we consider the following Bianchi equation: 



Since Vp, ip, p, a, (3 can be bounded in L™L™L°°(S) } it follows from Proposition [7] and Gron- 
wall's inequality that (3 can be estimated in L°° (S) , as desired. □ 

5.10.3. Smoothness of Spacetime away from H_ u . We now show that the spacetime is smooth 
away from the null hypersurface u = u s . In the region u<u s before the impulse, this follows 
from standard theory of local-wellposedness. We focus on the region after the impulse, where 
u>u s . In the following Proposition, we prove uniform estimates for the V4 and V derivatives 
of cifoo in the region u > u s . The bounds for all derivatives of the curvature components follow 
from a combination of the estimates on a^. Theorem [9] and the Bianchi equations. 

Proposition 77. Suppose the data for the impulsive gravitational wave satisfy the estimates 
as given in Proposition^ for some I > 2, K > 0. Then 

Y Y limsup HV^aoollz^^j) < C' IK . 



88 JONATHAN LUK AND IGOR RODNIANSKI 

Proof. By Proposition [33j 

Yl Yl supHV^V'ttnll^ir^+a-v)) < C' IK . (69) 



Fix 5 > 0. Take n large enough such that 2 11 < 8. Then (69) implies that 

S ^V\Nt^ ta n\\LHH u (u s +5,e))<C' ItK . 
i<I-l k<mm{K,[^=±\} % 

The conclusion follows from the fact that this bound is uniform in 5. □ 

As a consequence, we show 
Proposition 78. The limiting spacetime is smooth away from the the null hypersurface H_ u . 



Remark 9. In Proposition 71 , we have only used the fact that for u > u s , the initial data 
set is smooth. Therefore, suppose we have a initial data set satisfying the assumptions of 
Theorem [3J with the additional assumption that it is smooth for u>u for some u, we can 
also prove that the spacetime is smooth in < u < t, u < u < e. 

5.10.4. Decomposition of into Singular and Regular Parts. Finally, we show that has 
a delta singularity on the incoming null hypersurface H_ u . 

Proposition 79. can be decomposed as 

«oo = 5(u s )a s + a r , 

where S(u s ) is the scalar delta function supported on the null hypersurface H_ u , a s = 
a s (u,$) ^ belongs to L™L°°(S) and a r belongs to L™L™L°°(S). 

Proof. Define 

cts^^) := lim VT l x(u, u, i3) — lim Q,~ x(u, u, 

and 

ct r := - 5(u s )a s . 

We now show that a s and a r have the desired property. By Theorem [6j a s belongs to 
L^°L 00 (S'). That a s ^ follows from the fact that \ has a jump discontinuity across u = u s , 



which is proved in Proposition 75 



It remains to show that a r belongs to L°° (S) . To show this, we consider the measure 

of the half open interval [0,w) using the measure a r (u, 

(o r (u,i?))([0,u)) 

= (fi- 1 x)(M,M,^) - lim (IT 1 *) + lim (fT^Xw, u, ti) - (fT^Xw, u = 0, i?) 



Jo 

f- d f- d 

lim / — (f2 _1 x)(M, u, $)du + lim / — [pr x){ u i u, $)du 
Jo ou u-^uf Ju cm 



u-m s Jo 

+ / (toXX)(u,u,$)du. 
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By Proposition 
H u . Thus 
function 



77 



A. 

du 1 



(Q 1 x){u, u, $) is a bounded function away from the the hypersurface 
a r (uPd))([0,u)) can be expressed as an integral over [0,u) of an L°° (S) 

as desired. □ 
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